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SUMMARY 


This report presents the results of a program of analysis and test 
to investigate the effects of temperature, dissipation level and 
geometry on the dynamic behavior of elastomer elements. The report 
presents, first, a review of force displacement relationships in 
elastomer elements and the effects of frequency, geometry and 
temperature upon these relationships. Based on this review, methods 
of reducing stiffness and damping data for shear and compression test 
elements to material properties (storage and loss moduli) and 
empirical geometry factors are developed and tested using previously 
generated experimental data, A prediction method which accounts for 
large amplitudes of deformation is developed on the assumption that 
their effect is to increase temperature through the elastomers, there- 
by modifying the local material properties. Various simple methods 
of predicting the radial stiffness of ring cartridge elements are 
developed and compared. 


The test method for all dynamic stiffness and damping measurements 
involves resonating a mass, supported on the test element, by exciting 
the base to which the test element is mounted. The specimens tested 
were: Two rectangular shear specimens; one cylindrical compression 

specimen; and two cartridge elements of different length. Tests were 
performed at temperatures between 5 and 93 C and at dissipation levels 
between 0.022 watts/cm^ and 4.4 watts/cm^ material. 


Material properties were determined from the shear specimen tests as 
a function of frequency and temperature. Using these material 
properties, numerical predictions of stiffness and damping for 
cartridge and compression specimens were made and compared with corre- 
sponding measurements at different temperatures, with encouraging 
results. 


Successful comparisons of high dissipation predictions and measure- 
ments were made for the compression specimens but it is shown that, 
for the shear specimen, strain is an important influential parameter, 
even when the temperature rise is small. 

The method of reduced variables (which allows frequency dependent 
data at different temperatures to be collapsed to a single curve in 
reduced variables) is tested out on the material property results. 

It is shown that the storage modulus follows the reduced variable law 
satisfactorily whereas the loss modulus does not. 
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INTRODUCTION 


The use of support damping as a means to control rotating machinery 
vibrations is seeing increasingly wide application in advanced turbo- 
machinery. It will also have an important role to play in advanced, 
flexible, power transmission shafting. Presently the most common 
type of damper in these applications is the squeeze-film damper in 
parallel with some type of mechanical flexure. 

Elastomer dampers are an attractive alternative to the squeeze film 
for rotating machinery and other applications because of their sim- 
plicity; their inherent combination of stiffness and damping; their 
compactness; and their lack of need for seals or oil supply. In the 
form of 0-rings or cartridges they are being considered for low-cost 
engine applications and for helicopter transmission shafting. 

Two of the factors which resist the growth in application of elastomer 
dampers are the limited availability of design-oriented data on their 
dynamic behavior and limited quantification of the problems to be 
encountered in their application. Dynamic testing under controlled 
conditions, coupled with practical interpretation of the test results 
is the means to fill this need. Thereby, the influence of important 
geometrical, environmental and chemical design parameters can be 
determined . 

A program for testing and analysis of elastomer materials has been 
in progress at MTI since 1971. The objectives of the program are to 
develop elastomer dynamics technology and to reduce the technology 
to a form in which it can be readily applied in design. 

Experimental and analytical methods are employed in meeting these 
objectives. Firstly, from a limited number of tests, general prop- 
erties are sought which characterize an elastomer material. Secondly, 
prediction methods are sought which employ these material properties, 
in combination with appropriate geometry effects, to yield values 
for stiffness and damping of particular elastomer configurations. 

This report is the third in a series. Preceding reports (refs. 1 
and 2) have described the development and refinement of a powerful 
test technique, the Base Excitation Resonant Mass (BERM) method, 
and its initial application to the development of elastomer element 
properties in a controlled temperature environment under conditions 
of low dissipation. The BERM method most conveniently supplies data 
in the 100 Hz to 1000 Hz frequency range - into which the majority 
of rotating machinery vibration problems fall. Important features 
of the method are its controlled temperature environment; its ability 
to apply controlled high or low levels of dissipation at frequencies 
up to 1000 Hz or higher; its exploitation of mechanical resonance to 
achieve phase angles clustered around 90 degrees; and its determination 
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of stiffness and damping from accelerometer signals only, each 
processed by identical electronics, to minimize spurious phase 
errors in the data . 

The work presented herein provides a much expanded body of test 
data under a range of frequencies, temperatures and dissipation 
levels; development of interpretive and predictive methods; appli- 
cation of these methods; and verification of their ability to 
predict dynamic behavior. As a framework against which this work 
can be viewed, Figure 1 is presented. 

Figure 1 is a statement of the possible interactions between tests, 
data reduction, analysis, performance prediction and verification. 

As shown, the starting point for practical information is testing, 
and the BERM test method provides the main source of test data 
within this program. The output of the tests is raw data (ampli- 
tudes, phase angles, etc) which is converted by automated data 
analysis methods to specimen properties - that is, stiffness and 
damping of the elements being tested, as a function of frequency. 

To replace the discrete data points resulting from tests with a 
continuous relationship which can be handled mathematically, an 
appropriate frequency correlation, preferably of minimum complexity, 
is performed. The output of this correlation provides a compact 
definition of the dynamic behavior of the elements being tested 
which, together with similar compact definitions for other elements, 
yields a specimen catalog; that is, a list of stiffness and damping . 
properties of a series of different elements under different con- 
ditions of temperature and power level. 

To provide more universally applicable information, the specimen 
information is processed to yield material properties (which relate 
stress to strain) as a function of frequency. Where available, 
other sources of test data may also be employed. The output of this 
step is a material properties catalog, describing one material or 
several materials as a function of frequency and temperature. 

To apply elastomer material properties (moduli) in the design pre- 
diction of elastomer element properties, it is necessary to trans- 
late stress-strain relationships back into force - displacement 
relationships. For each elastomer configuration, there exists a 
ratio between the stiffness and damping at a particular frequency 
and temperature and the moduli of the material at the same frequency 
and temperature. These ratios, or geometry factors, will, in 
general, be a function of frequency, temperature, amplitude and 
geometry . 

As shown in Figure 1, sources of these ratios are solutions to the 
continuum equations for the elastomer element, in some cases com- 
bined with empirical coefficients obtained by tests. The spectrum 
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of solution complexity can run from a simple algebraic formula to 
the solution of complex continuum energy and viscoelasticity equa- 
tions in two or three dimensions. The body of ratios, or geometry 
factors, determined in one or several ways, is the geometry catalog. 

The material properties catalog and the geometry catalog are the two 
sources of information which embody the predictive model. They make 
it possible to predict stiffness and damping as a general function 
of material, geometry, frequency, amplitude and temperature. 

Progress along most of the paths of Figure 1 has been made under 
this program, as will be described in this report. 

Important questions to be addressed by this technology program are: 

”What level of complexity and sophistication is needed to 
generate an adequate geometry catalog for design purposes?** 

**What are the benefits, in improved accuracy, to be expected 
from advanced solutions relative to simple solution methods?** 

The simplest methods of all provide geometry factors which are in- 
dependent of frequency, and in two cases such predictive methods 
have been tested and their accuracy assessed against test data or 
more complete representations. These methods clearly lack rigor and 
there are discrepancies between their predictions and measured data. 

At the same time, for certain application of elastomers, their pre- 
cision may be adequate. In other applications, the development of, 
as yet, unavailable sophisticated prediction methods will be required. 
Given the level of component prediction accuracy to be expected, as 
indicated in this report, sensitivity studies for particular applica- 
tions will reveal the likely uncertainties in predicted system per- 
formance which will result. 

The following sections of the report present, first, a summary of 
results, conclusions, and recommendations, including the material 
properties and geometry catalog resulting from the tests and anal- 
ysis presented in the report. Then follows a review of force- 
displacement relationships for elastomers; the development of 
interpretive and predictive models; a description of test methods, 
including parameter ranges; a presentation of the test results; and 
a discussion of these results. 
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SUMMARY OF RESULTS 


1 . 


2 . 


3. 


4. 


5 . 


6 . 


7. 


8 . 


9. 


10 . 

11 . 


12. 


Methods of identifying generalized material properties and 
their relationship to dynamic behavior of elastomer elements 
have been developed, applied and tested. 


Dynamic test data for several geometries has been generated, 
for dissipation levels from 0.022 to 4,4 watts/cm'^? for tem- 
peratures from 5 to 93 C, and for frequencies from 70 to 
1500 Hz. 


The geometries considered are: 

• shear loaded rectangular elastomer strips, 

• compression loaded right solid cylinders, and 

• radially loaded cylindrical cartridges. 

The material is a broad temperature range (BTR) elastomer- 
polybutadiene . 

Material property values (storage and loss moduli) have been 
generated, as a function of frequency and temperature, from 
the shear elements. 

Geometry factor values (analogous to the traditional shape 
factor coefficient) have been generated for compression loaded 
right circular cylinders as a function of frequency. 

These geometry factors have been shown to provide a self- 
consistent representation of data generated under previously 
reported work. 

These geometry factors have provided reasonably effective 
predictions of stiffness and damping for a range of temperatures 
other than those for which they were developed. 

These empirical geometry factors, which depend on frequency, 
have been compared, in terms of their ability to predict 
measured stiffness and damping, with use of a frequency in- 
dependent static shape factor coefficient . 

Cartridge test specimens and an appropriate test fixture have 
been designed and fabricated. 

Predictions of cartridge stiffness and damping hafve been made 
as a function of frequency, using two different algebraic 
expressions. The simplest of these is based on representation 
of the cartridges as a series of noninteracting beams and 
column elements. 

These two simple methods predict results which differ by a 
factor of two but provide a band within which the great 
majority of measured values lie. 
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13. A thermoviscoelastic prediction method indicates some similar 
trends for stiffness as observed under conditions of high dis- 
sipation; that is, decreasing stiffness with increasing dis- 
sipation levels. This statement applies particulary to compres- 
sion test specimens. 

14. Although the conductivity is unknown for this material, an effective 
thermal conductivity value has been generated which gives a 

good, quantitative match between measured and predicted values of 
high dissipation stiffness and damping for compression loaded cylinders. 

15. Use of the same effective conductivity value for shear elements 
gives predicted stiffness and damping values which are higher 
than observed, under high dissipation conditions. 

16. Measured stiffness and damping values were correlated against 
strain and against dissipation. More consistent dependence 
on strain than on dissipation was observed, 

17. A ‘'strain softening" effect is apparent in the shear specimen 
results . 

18. The method of reduced variables was effective in reducing 
storage modulus data for different temperatures to a single 
curve of reduced modulus vs reduced frequency, but was not as 
effective for loss modulus. Strain dependence of the results 
is shown to be a likely cause of this apparent inconsistency 
for loss modulus . 

19. The material property catalog for the tested elastomer, based 
on the results presented in this report is as follows: 


Temper- 

ature 

1 Storage Modulus I 

Loss Modulus 

32 C 
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o 

r— 1 
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0) 
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G^(o)) 
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03 
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10^ 

-0.0637 

0) 
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G'(o)) 
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10^ 

0.3037 

0) 

G^1[o3) 

= 2.423 

X 

10^ 

-0.0301 

U) 
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= 1.197 
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10^ 

0.2450 

0) 

G"(o3) 

= 1.064 

X 

10^ 

+0.2246 

03 


20. The geometry catalog for the three tested goemetries is as 
follows : 
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CONCLUSIONS AND RECOMMENDATIONS 


On the basis of the results presented in this report, the following 
conclusions are drawn. 

1. Storage and loss modulus values, obtained as a function of 
frequency from shear specimen tests, are an effective 
representation of the dynamic properties of a material. 

2. Empirical, frequency dependent coefficients provide a self- 
consistent representation of the variation of stiffness with 
shape factor for cylindrical compression elements. 

3. For new conditions, the frequency dependent empirical co- 
efficients provide similar predictive ability to the use of a 
constant, static, coefficient. 

4. Either method of prediction for the compression elements 
incurs errors of 10 to 30 percent (high) relative to test 
data generated under the present program. 

5. Predictive ability, based on simple static analyses, for 
radially loaded ring cartridges provides upper and lower limits 
for stiffness and damping within which the great majority of 
measured data lies. 

6. The effect of high dissipation is to reduce stiffness. 

7. According to the frequency and geometry damping can be re- 
duced or increased by high dissipation. 

8. Self-heating and strain-softening effects are the combined 
causes of stiffness and damping changes under high dissipa- 
tion . 

9. For all results presented, stiffness and damping depend more 
consistently on strain than on dissipation rate. 

10. Increasing strain up to 0.02 (double amplitude) generally in- 
creases the loss coefficient for elements of polybutadiene. 

For future investigations, the following recommendations are made. 

1. The understanding of the dynamic behavior of elastomer elements 
will be improved by obtaining further data at very low strain 
levels, and by relating it to the present data at high strains. 

2. Simple methods for including strain dependence of dynamic be- 
havior in elastomer design practice should be developed. 
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3* The thermal properties of the elastomer material under investi- 
gation should be determined. 

4. Material and geometrical property tests for a number of dif- 
ferent elastomer materials should be performed to determine the 
extent to which the presently reported results can be 
generalized . 

5. Predictive ability for radially loaded cartridges must be 
further evaluated in a rotating load environment. 

6. Because they are compact, readily available, elastomer 
elements, the dynamic properties of 0-rings should be investi- 
gated by analysis and by test, with an effort being made to 
relate the dynamic properties to the basic material properties 
of storage and loss moduli. 

7. Application studies, analytical then experimental, should be 
performed to determine the benefits and limitations of 
elastomer elements as vibration control elements, 

8. Successful short-term application tests should be followed by 
tests for life and survivability to determine the extent to 
which these problems will limit application. 

9. Application studies should be interpreted to reveal what level 
of predictive accuracy is required for elastomer elements. 

10. More sophisticated prediction methods should be developed as 
indicated by the application studies. 
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REVIEW OF MATERIAL PROPERTIES AND THE EFFECTS 
OF GEOMETRY, FREQUENCY, AMPLITUDE AND TEMPERATURE 


In this section the relationships between force, displacement, stress 
and strain in an elastomer specimen are reviewed, with particular 
emphasis on frequency domain behavior and the influence of geometry, 
temperature and amplitude. 

Relaxation and Hysteresis 

Elastomer deformation properties are dominated by the process of re- 
laxation. By this process, the deflection for a given load (strain 
for a given stress) increases as a function of time - tending towards 
an asymptotic limit after a period of sufficient length. For very 
short periods of time, the displacement under a particular load will 
be significantly lower than the sympototic limit. 

If a load cycle is applied repeatedly to an elastomer element, the 
process of relaxation causes hysteresis - that is, a load deformation 
curve which depends on the rate of change of load as well as on the 
load's magnitude - see Figure 2, in which a representative elliptical 
hysteresis loop is presented. Implicit in such a load deflection 
curve or hysteresis loop is a loss of energy per cycle equal to the 
area of the hysteresis loop. Expressed in vibration terminology, 
elements with such a load-deflection curve provide dissipation, or 
damping. As shown by Lazan (ref. 3), for an elliptical hysteresis 
loop, the material behavior is linear, and the damping may be expressed 
as 


D 

K (co) = — ~ (1) 


where D is the area of the hysteresis loop in energy units and Y^ is 
the sinfle amplitude. It is also possible to define stiffness in 
terms of this loop as the ratio between force at maximum delfection 
and the value of maximum deflection. 

K^(co) = F^/Y^ (2) 

where and Y are defined in Figure 2. In the following paragraphs 
the relationships between instantaneous force and displacement are 
developed . 

It is noted, at this point, that the linearity of elastomer be- 
havior has been investigated by Warnaka (ref. 4) for materials 
including neoprene, nitride rubber and BTR elastomer for strains 
up to 0.1. Using electronics to sense deviations from an el- 
liptical hysteresis loop, Warnaka concluded that linear dynamic 


11 



i FORCE, F 



I 


Load-Deflection Curve For an Elastomer Showing Hysteresis as a Result of 
The Relaxation Process 




behavior within a small discrepancy was exhibited by elastomer elements 
for even the highest strain at frequencies in the 100-1000 Hz range. 


Relationships Between Force And Displacement in The Frequency Domain 

A linear elastomer element undergoing a repeated deflection cycle 
defined by 


Y 




(3) 


gives rise to a 


cyclic force variation described by 


F 


i03t 

F e 


(4) 


where: cd 

Y 
F 


* 

Y 



is the frequency (rad/sec) 
is the instantaneous displacement (m) 
is the instantaneous force (N) 
is time 

are complex numbers defining the amplitude and phase 
of the displacement and force variations, (and it is 
understood that only the real parts of expressions 
(3) and (4) are used to establish values for Y and F) . 


The force and displacement may be related by 

* 

F = K Y 


where K 


is a complex number 


•k 

K 


+ iK^ 


(5) 

( 6 ) 


and is, generally, a function of frequency. The quantities Kj^ and 
K will be referred to as the stiffness and dmaping of the elastomer 
element. The and K_, which depends on frequency, temperature 
amd amplitude characterize the dynamic behavior of the element. The 
damping ability of an elastomer element is sometimes expressed in 
terms of a loss coefficient (nondimens ional) which relates damping 
to stiffness, or the energy dissipated percycle (D^) to the energy 
stored at maximum displacement (U^) 


" " (co) 2 n Ug 


(7) 


Relationships Between Stress And Strain 

Considering any elastomer clement, its geometry can, normally* be 
characterized by a certain area of the element (A) which is stressed 
and a certain dimension of the element, h, which is strained. 

Figure 3 shows two elastomer element geometries which have been 
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tested for dynamic stiffness and damping, and which illustrate the 
effects of geometry. For the shear element (Fig. 3A) , the stressed 
area. A, is its width times its length (b*x) and the strained dimen- 
sion is its thickness, h. For the compression element (Fig# 3B), 
the stressed area is the top area of the cylinder (n D^/4) and its 
strained dimension is the height, h. 

For the shear element stiffness and damping can be related to effec- 
tive moduli as follows. 


Stiffness : 

K, = G%. ^ 
1 ef f h 

(8) 

Damping 

K = g'' , ^ 
2 eff h 

(9) 


where sre shear moduli, commonly referred to as the 

effective storage and loss moduli, respectively. 

For the compression element stiffness and damping can be related to 


the effective moduli as follows 





Stiffness : 


= 3 

eff h 

(10) 

Damp ing : 

K 

2 

= 3 

g " - 

eff h 

(11) 


where the factor of 3 arises from the classical elasticity relation- 
ship E = 2(1 + v)G and the fact that elastomer Poisson's ratios are 
typically in the range 0,4945 to 0,4999 (Holownia, ref. 5). Any 
errors in applying this classical elasticity relationship to dynamic 
behavior of a viscoelastic material will be implicit in the state- 
ment that ^^^ff effective moduli. ^eff generally 

functions or the material, the frequency, the amplitude of loading, 
the temperature, the initial strain and the geometry . They are, 
therefore, component rather than material properties. Further dis- 
cussion will explore the relationships between effective and "true" 
moduli . 

The Effects of Geometry 

For static loading the effects of bending can be shown (ref. 6) to 
cause the following relationship between effective and true modulus 
for the shear element of Figure 3A. 


[l + h2/3x2 ] 


( 12 ) 


where X is the dimension of the sheared area in the direction parallel 
to the loading, h is the thickness, and the subscript, o, implies 
static loading. In Figure 4, the stiffness correction factor 
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(G^f^/G^) is plotted as a function of X/h for a rectangular shear 
element. Under static compression loading, the interaction of axial, 
radial and circumferential stresses was shown by Gent and Lindley 
(ref. 7), Hattori and Takei (ref. 8), and Payne (ref. 9) to give the 
following relationship between apparent, or effective static modulus 
and actual, or true, static modulus: 

KifK ' “o' 

where S is a shape factor defined as the ratio of the loaded area 
to the unloaded area (D/4h for a cylinder) and p is a coefficient, 
shown by Gent and Lindley (ref. 7) to be equal to 2, 

A more complete expression based on the work of Gent and Lindley 
(ref. 7), which attempts to account for bulk compressibility effects 
is : 


9 1 3 G' 

7p ^ N = (1 + (14) 

^ eff^o ® 

where B is the bulk modulus for the material. However, this ex- 
pression has been tested by Moghe and Neff (ref. 10) using classical 
elasticity solutions and shown also to be incomplete for large 
values of S in that an additional dependence on Poisson's ratio 
occurs. For values of S less than 10, however, neither the second 
term in Eq . (14) nor the inadequacies due to Poisson's ratio effects 
should be significant. 

In Figure 5, the static stiffness correction factor ((Ggf f 
as given by Eq . (13) with p =2, is plotted against the ratio of 
diameter to height, for a compression loaded cylinder 

The existence of such a square law dependence on shape factor has 
been verified experimentally by Payne (ref. 9). He showed that, for 
natural rubbers, the shape coefficient, p, was close to 2 but, for 
harder materials the coefficient p could fall by as much as 50 per- 
cent . 

Under low frequency dynamic loading (10 Hz) Payne also demonstrated 
that the effective storage modulus varied with shape factor according 
to a square law. 

The analytical or experimental quantification of the effects of 
geometry (shape) at higher frequencies appears to be limited. Nashif, 
Cannon and Jones (ref. 11) have demonstrated that shape has a pro- 
nounced influence on the effective storage modulus of cylindrical 
specimens, in the frequency range of 250 to 500 Hz, but have not 
generalized their results. 
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The Effects of Frequency 

As might be expected from a material whose behavior is governed by 
time dependent relaxation processes, the slope and area of the 
hysteresis loop, together with the stiffness, damping, and loss 
coefficient are functions of the frequency with which the load 
cycle is executed. It is generally observed that the dynamic 
stiffness is higher than the static value, and that stiffness tends 
to increase with increasing frequency. Shown on a scale which may 
cover many decades of frequency, the general form of elastomer be- 
havior variation has three regions (ref. 12): 

• Rubbery 

• Transition (Leathery) 

• Glassy 

These are illustrated in Figure 6. 

In the rubbery region, the behavior is, as implied ’^rubbery.*’ Elastomer 
elements are resilient and the damping is usually of the order of 10 
to 30 percent of the stiffness. The hysteresis loop is a rather 
narrow ellipse. As the frequency increases, the transition region 
is encountered. Both stiffness and damping increase more rapidly in 
this region but the damping moves ahead of the stiffness and the 
hysteresis ellipse becomes broader. 

For a limited frequency range, damping (K 2 ) becomes higher than 
stiffness (Kj[) and the loss coefficient becomes greater than 1. For 
further increases in frequency the stiffness increases more than the 
damping and in the glassy region both quantities level off at values 
giving a loss factor which again tends to be 0.3 or less. 

The rubbery and transition regions are those normally encountered in 
elastomers and offer the best opportunity for use in vibration con- 
trol . 

Representation of Frequency Dependence 

Observed variation of stiffness and damping with frequency gives 
rise to discrete data points. For mathematical purposes it is con- 
venient to seek, for an element, a continuous relationship between 
force and displacement (or stress and strain) which closely matches 
observed behavior. There are several possible approaches to this 
problem, with different levels of sophistication and convenience, 
as discussed below. 

1, Generalized Viscoelastic Models . An approach to representing 
the variation of stiffness and damping of an elastomer element with 
frequency is to replace the elastomer with a mechanical analog, con- 
sisting of a series of springs and viscous dashpots. The simplest 
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of these models are the Maxwell and Voigt elements (Fig. 7). Using 
these basic elements as building blocks, more complex systems of 
springs and dashpots may be developed, such as those shown in 
Figure 8 . 

The application of a triparameter model, as shown in Figure 9, to 
elastomer specimen test data generated by the BERM method was dem- 
onstrated by Chiang, Tessarzik and Badgley (ref. 1). The methodology 
to determine, for a selected general viscoelastic model, the element 
values to best fit a particular variation of stiffness and damping 
with frequency was developed and demonstrated by Gupta, Tessarzik 
and Cziglenyi (ref. 2). 

2, Generalized Force - Displacement Relationships . An alternative 
approach for representing the variation with frequency is to assume 
a general linear differential relationship between force and dis- 
placement : 


2 3 

T7 ^ dF . d F . d^F , 

a F + a. 77+3 — r + a — r + 

° ^ 2 dt^ ^ dt^ 


„ , dY ^ d^Y ^ 

= c Y + c — + c — - + . . . 
o 1 dt 2 .2 

at 


(15) 


For sinusoidal loading and displacement: 


So that Eq . (15) becomes: 
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A close connection between this generalized force displacement re- 
lationship and the generalized viscoelastic models exists since, for 
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a given order of expansion on left- and right-hand sides of Eq . (15), 
there exists a corresponding combination of springs and dashpots 
(ref. 12), which will give the same force-displacement relationship. 

For example the differential equation for the triparameter model 
(Fig. 9) is 

^ " W ^ <-<12 + “ll> S 

so that the nonzero coefficients of Eq . (15) for the triparameter 
model are 


a = 1 

o 


Cti/Kti 


= K 


T2 


( 21 ) 


Cl ( K .^2 + 

3. Direct Fit of Mathematical Expressions . A direct approach to 
the problem is to impose upon the data some correlation function and 
to seek those coefficients which give the best match between the 
function and the data. Polynomial, trigonometrical, or power law 
(Kf = A^(aj relationships are some of the options available for this 
approach. It was found by Gupta, et al (ref. 2), that the use of a 
power law relationship was quite adequate to represent the behavior 
of elements of polybutadiene in the frequency range 100 to 1000 Hz. 
For engineering purposes, a known relationship between stiffness, 
damping and frequency for elastomer elements of different geometries 
is all that is required to use these elements in design for control 
of machinery vibrations. In essence this approach may be considered 
as the use of a generalized Voigt element in which the storage and 
dissipative elements are a continuous function of frequency. 


Influence of temperature . Elastomer behavior is, in general, 
strongly influenced by temperature. Thus complete characterizat ion 
of elastomer elements must provide their dynamic behavior over the 
frequency and temperature ranges of interest. In considering the 
effects of temperature the following observations are often made. 


1. Temperature tends to be most influential where 
frequency is also strongly influential. 

2, Increasing temperature tends to produce the same 
results as decreasing frequency and vice versa. 


Starting from these observations, an elegant method of reducing the 
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dynamic behavior of elastomers (and other polymers), covering a 
range of frequencies and temperatures, to a single curve relating 
appropriate reduced properties, was developed by Ferry, Fitzgerald, 
Grandine and Williams (ref. 13). The method is sometimes referred 
to as the WLF law, due to a paper by Williams, Landell and Ferry 
(ref. 14), in which the mechanics of its implementation are described 
and demonstrated. The mathematical statement of this method is as 
follows 








^300), 6c 


S.86 (e - 0 ) 

c 

101.6 + (9 - 9 ^) 
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( 22 ) 
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(24) 


where 0 is a characteristic temperature given by 
c 

6 


= 9^ + (50° ± 5°) K 


(25) 


and 6 


is the transition temperature between the rubbery 

and glassy states of the material. 

is the observed modulus at temperature 9 and 

frequency o), 

is the reduced modulus, 

is the coefficient by which test frequency is 
multiplied to obtain the corresponding reduced 
frequency at temperature 0^ , 


By these relationships, knowledge of a modulus as a function of fre- 
quency at several temperatures may be interpreted as knowledge of 
reduced modulus over a wide frequency range at the characteristic 
temperature. This knowledge may be further interpreted to give the 
variation of modulus over a wide frequency range at any temperature 

between 9 and the maximum temperature tested, 
c 

The method has been demonstrated to be applicable to material prop- 
erties of storage and loss moduli or to element properties of stiff- 
ness and damping. The additional information needed to implement 
and use the methods, beyond the elastomer dynamic data is the value 

of 0 . 
c 


The transition temperature, 9t, may be measured independently or 9 
may be inferred from a sufficiently large volume of dynamic data. 
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Influence of large amplitude displacements . Observed properties 
of an elastomer element tend to be a function of amplitude of 
loading, as well as frequency and temperature. This is shown, for 
example, by references 15, 16 and 17. No specific and complete 
explanation of this mechanism has been produced. It may be 
hypothesized, however, that a major contribution is the increase 
in temperature of the element when a significant dissipation rate 
is imposed. As will be discussed, further, in relation to the test 
data a second, nonthermal, mechanism has also been shown to cause 
degradation in stiffness at high strains. 
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INTERPRETATION OF BASE EXCITATION RESONANT MASS TEST DATA 


The preceding review of elastomer material properties was presented 
as a framework against which the methods of dynamic testing described 
in this report and the methods of data reduction and interpretation 
may be viewed . 

The BERM method has been applied in the generation of dynamic prop- 
erties of three categories of elastomer element: shear elements; 

compression elements; and ring cartridge elements. The shear and 
compression elements have been illustrated in Figure 3. The ring 
cartridge element is shown in Figure 10. 

Table I summarizes the range of conditions covered by this data and 
the periods when it was generated. The 1973, Phase II, data was 
previously reported by Gupta, Tessarzik and Cziglenyi (ref. 2) in 
terms of stiffness and damping as a function of frequency. To 
obtain a continuous mathematical representation of the dependence 
of stiffness and damping on frequency, a power law relationship was 
fitted to the test data and the results of this data analysis are 
presented in Table II. 

Prior to the generation of the Phase III test data presented in this 
report, an effort was undertaken to develop methods of exploiting 
the information generated by the BERM test method, to yield material 
properties and geometry factors which can be applied in engineering 
design. The data of Table II was used to implement and test these 
methods . 

Generation of Material Properties 

In the previous section, it was shown that stiffness and damping 
values can be interpreted in terms of relationships between stress 
and strain. It was also shown that combined stress effects tan be 
instrumental in making the effective moduli for the shear and com- 
pression element different from uniaxial moduli. Table III sum- 
marizes the static stiffness correction factor 

would be expected on the basis of Eqs. (12) and (13) for the shear 
and compression elements tested under Phase II and III. 

In all but two of the elements tested the predicted effects of 
combined stress are greater than ten percent, l^ile this data is 
for static behavior of the elements, and no rigorous basis exists 
for applying it to dynamic behavior, it is assumed that the two 
shear elements for which the static stiffness correction factor is 
within two percent of unity are also likely to show small combined 
stress effects under dynamic loading conditions. All the other 
elements tested are Likely to show significant coinbinecl stress 
effects . 
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TABLE I 


SUMMARY OF DYNAMIC TEST CONDITIONS 
Material Tested - Polybutadiene 


Phase 

Da te 

Geometry 

Dimensions (cm) 



Nomina 1 
Dissipa t ion 

W/cm^ 

h 

X 

b 

.... 

d 

D. 

i 

D 

0 

L 

II 

L973 

Shear 

0.3175 

0.508 

4.877 





KBI 


0.022 




0.3175 

1.27 

4.877 





8 

32 C 

0.022 




0.3175 

2.54 

4.877 





4 

32 C 

0.022 



Compression 

0.3175 



1.27 




3 

32 C 

0.022 




0.635 



1.27 




10 

32 C 

0.022 




1.27 



1.27 




30 

32 C 

0.022 

III 



0.3175 

1.27 

4.877 





8 

32 C 

0.062 




0.3175 

1.27 

4 .877 





8 

50 C 

0.062 




0.3175 

1.27 

4.877 





8 

66 C 

0.062 




0.3175 

1.27 

4.877 





3 

93 C 

0.062 




0.3175 

2.54 

4.877 





4 

32 C 

0.033 




0.3175 

2.54 

4.877 





4 

32 C 

0. 176 




0.3175 

2.54 

4.877 





4 

32 C 

0.352 




0.3175 

2.54 

4.877 





4 

32 C 

0.704 




0.3175 

2.54 

4.877 





4 

32 C 

1.76 




0.3175 

2 .54 

4.877 





4 

32 C 

4.40 


1974 

Compression 

0.635 



1.27 




10 

5 C 

0.062 




0.635 



1.27 




10 

66 C 

0.062 




0.635 



1.27 




10 

93 C 

0.062 




0.635 



1.27 




10 

32 C 

0.44 




0.635 



1.27 




10 

32 C 

1.] 




0.635 



1.27 




10 

32 C 

2 . 2 


1975 

Cartridge 





1.905 

2 .358 

0.476 

■ ■- 
4 

32 C 

0.022 








1.905 

2 .358 

1 0.476 

4 

43 C 

0.022 








1.905 

2 .858 

0.476 

4 

66 C 

0.022 








1.905 

2 .858 

0.476 

4 

80 C 

0.022 








1.905 

2 .858 

0.476 

4 

32 C 

0.11 








1-905 

2 .358 

0.476 

4 

32 C 

0.22 








1.905 

2 .858 

0.476 

4 

32 C 

0.88 








1.905 

2 .858 

0.953 

4 

32 C 

0.044 








1 .905 

2.858 

0.953 

4 

43 C 

0.044 








1 .905 

2.858 

0 . 953 

4 

66 C 

0.044 








1.905 

2.858 

0.953 

4 

80 C 

0.04 4 








1.905 

2.353 

0.953 

4 

32 C 

O.H 








1 .905 

2 .353 

0.953 

4 

32 C 

0.7 








1-905 

2 .353 

0.953 

4 

32 C 

1 .01 








1.905 

2 .358 

0.953 

4 

32 C 



I 

I 


I 




















TABLE II 


SUMMARY OF ELASTOMER TEST DATA FROM PHASE II TESTS 


MATERIAL • Polybutadiene 

TEMPERATURE « 32 deg-C 


CuMPLEX STIFFNESS K • Kl ♦ 1K2 ib/ln. 

CORRELATIUN FORM Y * A*OMEGA**8 m rad/sec) 


SPECIMEN 


*********«-**** Kl ****«*****^*ikiH^ 

*#********«#** K2 «*********#'*«* 


AMPL I TUDE *ir***itit**ik 

DETAILS 


MEAN 

A 

B 

MEAN 

A 

B 

MEAN 

A B 



Ib/ln 



ib/ln 



In 

0.50IN-30/C/2.5 

/.54IE^04 

2.224E^04 

I.589E-0I 

8.347E+03 

1 .873E+04 

-1 .05IE-0I 

1 .055E-03 

3.854E-OI -7.677E-OI 

0.50IN-30/C/5.0 

7.809E^04 

2.545E+04 

I.458E-0I 

8.657E^03 

I.507E+04 

-7.2I2E-02 

I.042E-03 

3.401 E-01 -7.529E-01 

0.25IN-I0/C/2.5 

d.042E^04 

2.87IE^04 

I.325E-0I 

9.266E4>03 

2 .809E+04 

-I.427E-0I 

3.806E-04 

2.040E-01 -8.082E-01 

0.25IN-I0/C/5.0 

e.026E^04 

3.577E+04 

1 .044E-0I 

9.447E+03 

2.753E+04 

-I.382E-OI 

3.924E-04 

1.764E-01 -7.889E-01 

0.I25IN-3/C/2.5 

7.662E4-04 

7.3I2E+04 

6.368E-03 

1 . 1 90E^04 

1 .378E+04 

-2.004E-02 

2.501 E-04 

2.407E-02 -6.2I6E-0I 

0. I25IN-3/C/5.0 

0.3 926^04 

7.007E+04 

3.227E-02 

1 .295E+04 

1 .854E+04 

-4.859E-02 

2.250E-04 

3.547E-02 -6.854E-0I 

0.5IN-8 

/S/0. 

I.798E+05 

5.330E^04 

1 .548E-0I 

1 . 720E^04 

3.457E+04 

-8.886E-02 

3. I86E-04 

2.478E-0I -8.473E-OI 

0.5 IN-6 

/5/2.S 

1 . 666E^05 

5.003E^04 

I.530E-0I 

1.737E>04 

3.I38E4-04 

-7.521 E-02 

3.498E-04 

4.985E-01 -9.232E-01 

0.5IN-8 

/S/5.0 

1 . 598E^05 

5.361 E*04 

1 .39 IE-01 

1 .653E+04 

2.678E+04 

-6.142E-02 

3.462E-04 

2.I62E-OI -8.I98E-0I 

1 .OIN-4 

/S/0. 

1 .453E405 

2.569E^04 

2.199E-0I 

l.624E^04 

2.800E4*04 

-6.9I2E-02 

3.497E-04 

2.844E-0I -a.503E-0l 

1 .OIN-4 

/S/2.5 

1.35IE-K)5 

2.550E4-04 

2.I35E-0I 

I. 63 I E^04 

2.388E^04 

-4.882E-02 

3.799E-04 

2.502E-0I -8.309E-0I 

1.0 IN-4 

/S/5.0 

\ . 307E+05 

2.263E^04 

2.238E-01 

I.563E+04 

I.923E+04 

-2.647E-02 

3.820E-04 

2.635E-0I -8.339E-0I 

0.2 IN-20 

/S/0. 

l.625E<^05 

3.708E+04 

I.881E-0I 

l.689E^04 

5.700B404 

-1 .548E-0I 

3.395E-04 

2.l68E-Oi -8.220E-0I 

0.2 IN-20 

/S/2.5 

1 . 549E^05 

3.336E+04 

I.936E-0I 

l.68IE^04 

5.867E^04 

-1.576E-0I 

3.398E-04 

2.630E-0I -0.386E-OI 

0.2IN-20 

/S/5.0 

I.452E405 

3.689E+04 

1.743E-01 

I .600E+04 

4.833E^04 

-I.406E-01 

3.683E-04 

2.240E-0I -8.I56B-01 





CUMPLEX 

STIFFNESS K 

- Kl ♦ 1K2 

New/m 






CORRELATIUN FORM Y 

« A*0MEGA**B 



SPECIMEN 

************** Kl *****★*****«*« 

Mr****** ****** K2 ************** 

************ 

AMPLITUDE ********** 

DETAILS 

MEAN 

A 

B 

MEAN 

A 

B 

MEAN 

A B 


New/m 



Ne»</m 



m 


0.50IN-30/C/2.5 

1 . 32 1 E+07 

3.894E+06 

I.589E-0I 

I.462E*06 

3.280E+06 

-1.05IE-0I 

2.679E-05 

9.790E-03 -7.677E-01 

0.501N-30/C/5.0 

1 . 368E+07 

4.457E+06 

I.458E-0I 

1.5I6E*06 

2.640E^06 

-7.212E-02 

2.646E-05 

8.639E-03 -7.529E-0I 

0.25IN-I0/C/2.5 

1.408E^07 

5.027E+06 

1 .325E-OI 

I.623E+06 

4.920E+06 

-I.427E-01 

9.668E-06 

5.I80E-03 -8.082E-0I 

0.25IN-I0/C/5.0 

1.406E^07 

6.264E^06 

1 .044E-01 

».654E*06 

4.822E^06 

-1 .382E-OI 

9.968E-06 

4.480E-03 -7.889E-01 

O.I25IN-3/C/2.5 

1 .342E^07 

I .28IE>07 

6.370E-03 

2.084E+06 

2.4I4E*06 

-2.004E-02 

6.354E-06 

6, 1 14E-04 -6-2I6E-0I 

0. I25IN-3/C/5.0 

l.557E^07 

1 .227E^07 

3.227E-02 

2.268E*06 

3.247E+06 

-4.059E-O2 

5.715E-06 

9.010E-04 -6.854E-01 

0.5JN-8 /S/0. 

3.l49E-f07 

9.335E-*-06 

1.54BE-01 

3.012E^06 

6.054E<-06 

-d.886E-02 

8.092E-06 

6.294E-03 -8.473E-0I 

0.5IN-8 /S/2.5 

2.918E^07 

8.762E+06 

1.530E-0I 

3.042E^06 

5.497E406 

-7.52IE-02 

8.805E-O6 

I.266E-02 -9.232E-01 

0.5IN-8 /S/5.0 

2.799E+07 

9.389E^06 

I.391E-0I 

2.895E^06 

4.609E^O6 

-6.I41E-02 

0.794E-O6 

5.491 E-03 -0. )98E-01 

I .OIN-4 /S/0. 

2.545E+07 

4.499E+06 

2. 199E-01 

2.844E^06 

4.903E*06 

-6.9I2E-02 

0.08 IE-06 

7.223B-03 -8.503E-0I 

1 .OIN-4 /S/2.5 

2.366E^07 

4.465E+06 

2.135E-0I 

2.057E*O6 

4.183E*06 

-4.883E-02 

9.648E-06 

6.356E-03 -8.309E-0I 

1 .0 IN-4 /S/5.0 

2.289E^07 

3.963E+06 

2.238E-0I 

2.737E-»-06 

3.369E^06 

-2.647E-02 

9.703E-06 

6.692E-03 -S.339E-01 

0.2IN-20 /S/0. 

2.S46E+07 

6.494E^06 

1 .88IE-0I 

2.959E*06 

9.984E*06 

-I.548E-OI 

8.624E-06 

5.506B-03 -8.220E-0I 

0.2IN-20 /S/2.5 

2.713E^07 

5.843E*06 

1.936E-01 

2,944£*06 

1.028E^07 

-I.576E-0I 

8.631 E-06 

6.680E-03 -8.386E-01 

0.2IN-20 /S/5.0 

2.542E+07 

6.461 E^06 

1 .743b-0I 

2.803b>06 

8.464E*06 

-1 .406E-0I 

9.355E-06 

5.691 E-03 -8. 156E-0I 

N.H. Frequoncy 

UMkGA In th# 

corr«intlor 

1 function 

H In rad/sec. 





A sample consists of several specimens of elastomer arranged in a prescribed 
configuration. 

Specimen details ' code given in the table is derived as follows^ 

(L) Compressit>n Sample: (Specimen tliickness) IN - No. of Specimens/C/ 

(7o Preload). 

The diameter of all compression specimens is 1/2 in. 

(2) Shear Sample: (Specimen length) IN - No. of Specimens/S/ 

(7o Preload) 

The specimen thickness and width normal to 

direction of deformation for all shear specimens 
are respectively 1/8 in and 1.92 in. 
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TABLE III. ELASTOMER TEST ELEMENTS - EFFECTS 
OF COMBINED STRESSES 


Element Type 

Length 

(cm) 

No. 

Diameter 

(d,cm) 

Strained 

Dimension 

(h,cm) 

Width 
(b , cm) 

^^eff^o 

G 

0 

Shear 

2.54 

B 


0.3175 

4.877 

0.995 

Shear 

— 

1.27 

8 




0.980 

Shear 

0.508 

20 




0.884 

Compression 


30 

1.27 

1.27 


1.125 

Compression 


10 

1.27 

0.635 


1.5 

Compression 


3 

1.27 

0.3175 


3.0 


This assumption leads to the following approach to interpretation of 
the test data. 


1. From dynamic stiffness and damping data for the shear 
elements of 1.27 cm and 2.54 cm length (parallel to 
loading direction) the material properties of storage 
and loss modulus should be inferred - neglecting what- 
ever small effects of geometry occur. 


From these elements, values of storage and loss modulus 
at a particular frequency are calculated as follows: 


G^(o)) 


N A 


(26) 


K_ (cu) h 

g"(cd) = (27) 

N A 

With values of stiffness and damping available for more 
than one specimen at a particular frequency, the values of 
G‘’(o)), G‘’^(o)) are calculated which minimize the rms deviation 
between stiffness and damping predicted using the G (u>) , 

G (oj) and the measured stiffness value for the element. 
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2* From test data for the compression elements of different 
shape should be extracted empirical coefficients which 
account for the effects of geometry. 

With material property data available from shear specimen tests, 
values of shape factor coefficient are determined which minimize 
the deviation between predictions and measurements of compres- 
sion specimen stiffness and damping. Separate coefficient values, 
P and are determined for stiffness and damping so that 

2 

K^(oo) = g'(co) + 3' S^] (28) 

¥^(ay) = g"(co) + 3" S^j (29) 

Implementation of Methods of Data Interpretation 

Using the Phase II shear element data for three preloads and two 
length values (six sets in all) rms best fit values of and G^' 
were determined over the frequency range 100 to 1000 Hz and the 
resultant set of discrete G^, G^' frequency points were then refitted 
to power law expressions. As a result: 


G'(cd) 

= 4.136 X 


(30) 

G"(03) 

= 3.369 X 

10^ ^-0.0793 

(31) 


Using the Phase II compression element data for two preloads and 
three different heights, independent, best fit values of p^, p^^ were 
determined at discrete frequency values between 100 and 1000 Hz and 
their frequency dependence was then expressed as a power law relation- 
ship so that 

Ggff(03) = G'(oi) [l + 12.33 (^) ] (32) 

G'Jff(cu) = g"((d) [l + 1.726 03*^-02992 ^^2 ^ 

Verification of Methods For Determining Material Properties And 
Geometry Factors 

Based upon the preceding results, individual shear elements of poly- 
butadiene would be expected to exhibit the following values of stiff- 
ness and damping (neglecting bending effects) : 
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(cd) 


(34) 


= 4 . 136 X 


,_6 0,184 b X 

10 OD — — 

h 


^2 


1 in6 -0.0793 b X 

3.369 X 10 CD 

h 


(35) 


and compression elements of polybutadiene would be expected to exhibit 
the following values of stiffness and damping: 


V / \ / lA^ 0.184 n li.no-:. -0.2896 fD 

Kj^ (to) = 4.136 X 10 CO liir’ 12.33 o j J 

V f .\ o O^IQ lA^ -0.0793 ttD^ r, ^ 0.02992 / D 

K 2 (co) = 3.369 X 10 CO j^l + 1.726 co 


(36) 


(37) 


To determine the "goodness of fit" given by the above expressions 
with the measured data from which they were derived, predicted stiff- 
ness and damping of shear and compression specimens are presented in 
Figures 11 through 22. Three different frequencies, 100, 300 and 
1000 Hz have been selected and at each of these, stiffness and damp- 
ing has been plotted against a geometrical parameter representative 
of the difference among the specimens of a particular class. For 
the shear specimens, this parameter is the ratio of sheared length 
to thickness (X/xh), for which data over a range from 0.2 to 1.0 is 
available. For the compression specimens, the independent variable 
is the ratio of diameter to height (D/h), for which data over a range 
from 1 to 4 is available. 


The consistency between the predictions and the test data is generally 
reasonable, although scatter in the test data leads to individual 
discrepancies of up to 20 percent of the predicted value. The mean 
discrepancy is about ten percent of the predicted value. 

The results presented in Figures 11 through 22 are taken as providing a 
satisfactory check of the methods presented for modeling and data 
interpretation. The same approach to development of material property 
data from shear specimen tests will be applied to the Phase III data 
to be presented subsequently. This data covers a much wider temper- 
ature range and also reflects improvements in the test method which 
evolved during the Phase III tests. 

Alternative Correlations For Shape 

The work of Nashif, Cannon and Jones (ref. 10) provides data on the 
stiffness of cylinders in compression. Dynamic stiffness values for 
a range of diameter to height ratios between 0.25 and 4.85 are given 
at 250 and 500 Hz. The material tested is paracryl (a Buna Nitrile 
elastomer manufactured by UniKoyal). Just as for the polybutadiene 
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Fig. 18 phase II Correlation. Single Compression Element. 
Diameter 0.5 In., 300 Hz Stiffness 
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Fig. 19 Phase II Correlation For Shape - Single Compression 
Diameter 0.5 In. , 1000 Hz 










Phase II Correlation for Shape - Single Compression 
Element, Diameter 0.5 In., 300 Hz Damping 
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Correlation with Work of Nashif, Cannon and Jones (Ref. 
Dynamic Stiffness of Cylindrical Elastomer Test Specimen 



stiffness data, this set of data has been used to calculate shape 
factor coefficient values for the correlation: 

G"A - 2, 

Stiffness: + p j ( 38 ) 

and the results for 250 and 500 Hz are = 1.484 and 1.485, respec- 
tively. The corresponding results for polybutadiene are 1.45 and 
1.29. Comparisons between these correlations and the test data of 
reference 11 are plotted in Figures 23 and 24 in terms of the ratio 
of actual stiffness to stiffness predicted for uniaxial stress. 

This independent confirmation of the ability of a square law correla- 
tion between stiffness correction factor and diameter to height ratio 
to match test data adds support to the use of such an empirical 
approach for consolidating the present compression data. 

Use of a Constant Shape Factor Coefficient 

As discussed previously in this report, the most commonly used value 
of static shape factor coefficient (3), for cylindrical compression 
elements is two. A very straightforward method of predicting dynamic 
stiffness and damping could be based on this coefficient value, as 
follows 


2 - 2 

Kj (») . 3G'(<») [l + 2 (^ ] (39) 

2 - 2 

Kj (OJ) . 3G"(») [l + 2 © ] (40) 

No analytical justification exists for this application of static 
coefficients but, since it has the appeal of simplicity and generality, 
the likely errors in the approach have been tested. Superimposed on 
Figures 17 through 22 are the above expressions. At high frequencies 
the predictions of these expressions for stiffness are 50 percent high 
relative to the more representative lines based on frequency-dependent 
coefficients. At lower frequencies the discrepancy for stiffness is less 
and the discrepancy for damping is closer to ten percent under all 
conditions . 

Effects of Temperature And High Dissipation 

Since the Phase II data did not include temperature as a varied 
parameter it cannot be used to evaluate any means of representing the 
effects of temperature. However, in anticipation of the Phase III 
data lor high temperature and high dissipation, a model intended to 
tie the effects of high dissipation to the change in properties as a 
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result of temperature increases, generated by energy dissipation in the 
elastomer, has been developed. This model simply needs to have avail-* 
able values of G and G as a function of frequency and temperature 

G"(o3,T) (41) 

G'"^(cd,T) (42) 

The hypothesis of this model is that, when large amplitude vibrations 
are imposed upon an elastomer element, they give rise to dissipation 
of energy throughout the element. As a result there is a tendency 
for the elastomer to heat up. The temperature distributions in the 
elastomer are governed by the equations of thermal conduction and 
appropriate boundary conditions. As a result of the distribution of 
temperature in the elastomer, and the dependence of elastomer moduli 
on temperature there is a distribution of material properties through 
the elastomer. Stress equilibrium within the elastomer, accounting 
for this distribution of properties, governs the forces generated by 
the elastomer, and its apparent stiffness and damping. 

The implementation of such a model in two or three dimensions could 
be very complicated. However, to typify the results given by such a 
model, and to identify the major items of information needed, while 
minimizing the complexity, a one-dimensional development has been 
undertaken. 


An important, valid, assumption is that the time scale associated 
with the vibrational phenomena (less than 0.01 second) is short 
compared to the time scale associated with changes in t he thermal 
condition of the material (h^ p C /k, which is of the order of 40 
seconds and higher). ^ 


One-Dimensional Model For Large Dissipation Effects 

This model is directed at either a shear or a compression element in 
which uniaxial stress prevails. The temperature distribution is one- 
dimensional and the flow of heat is normal to the stressed area. The 
equations to be solved are: 


1. The energy equation, in terms of cycle average temperature, T 



P c 


p 3t 



(43) 


2. The dynamic stress-equilibrium equation for uniaxial stress 



(44) 


where U is the amplitude of displacement relative to the base 
of the elastomer element. 
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3. The dependence of properties on frequency and temperature 

G* = G* (co,T) (45) 

4, Boundary conditions 


T T^atZ-0,h (46) 

U = 0 at Z = 0 (47) 

5. For imposed displacement at the top of the element 

u = Y* at Z = h (48) 


6. For base excitation, the equation of motion for the supported 
mass 




Y -- Mco U = 


2 * 
Y 

o 


(49) 


where 


Y is the amplitude of displacement of top relative 
to base, 

Yq is base excitation amplitude (absolute), 

Y is a factor equal to one for shear elements, 
three for compression elements, 

M is the supported mass, 

CD is the base excitation frequency 
A is the stressed area, and 

is the temperature at base and top. 


The complete development of Equations (43) to (45) is presented in 
Appendix A. 


^Qj-l^tion of the Equations . As a result of the dependence of G* on 
temperature the system of equations is nonlinear. While steady-state 
results are of most interest, solution in the time domain has been 
selected as more convenient to implement than an iterative solution. 
The solution method is a modification of the second order Crank- 
Nicholson Method (ref. 17), in which nonlinear terms are evaluated by 
extrapolation forward one-half time-step. 

The energy equation is written to apply at the time t + At/2. 
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t + iit 


(T 


t+At 
At 


T^) 


{-v<^' HI f} 


= R 


t+At/2 


t+At/2 


(50) 


The left-hand side has a linear operator and consists of known 
quantities at time t and unknown quantities at time t+At. The 
right-hand side is unknown since it is a function of temperature at 
time t + At/2, but can be approximated by the linear extrapolation 
from known quantities: 

V4/2 “ 

The equation for U is solved at time t+At, once T^. . and the 
associated values of G are known. 


The time integration solution is continued until results for succes- 
sive time steps indicate that thermal equilibrium has been reached. 
The equilibrium distribution of temperature, material properties, 
and deflection will determine the behavior which will be observed 
once the initial thermal transients have died out. 


An interesting feature of the thermal transient behavior occurs when 
the base excitation boundary conditions apply and the frequency of 
excitation and supported mass are such that the element causes sub- 
resonant vibrations when the element is cold but vibrates above 
resonance when the element gets hot due to self-heating. In this 
case there is a wide swing in phase angle during the thermal tran- 
sient and an inability to hold subresonant conditions for levels of 
dissipation higher than a certain limit. Figure 25 illustrates the 
behavior to be expected for a range of different input amplitudes 
and this phenonmenon is discussed further in the Results and Discussion 
Sections of this report. Observed stiffness and damping of the ele- 
ment at any time during the time integration, including equilibrium, 
are calculated as 


Kj - p I» {g* (53) 
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TIME - SECONDS 


Fij^. 23 High Dissipation Predictions Influence of Input 
Amplitude. 1.316 Kg Mass, 380 Hz, 32 C Ambient, 
Compression Specimen, 10 Elements, 1.27 cm 
Diameter, 0.635 cm High 
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It is important to note that a number of material properties not in- 
cluded in the dynamic tests are needed to implement this analysis; in 
psfticular , density, specific heat and conductivity. 

In addition, if this analysis is to be applied with meaning to com- 
pression elements in which the stiffness correction factor is of order 
1.5 and above, some account must be taken of the stiffness correction 
factor. In the present case, the significant assumption has been 
made that the stiffness correction factors follow the functional 
dependence on frequency and geometry determined empirically for cylin- 
^^ical elements and that the increase in effective modulus should take 
place on a local basis, and be reflected in the equilibrium equation: 


A. n* M" 
&Z L eff dZj 


0 


(54) 
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PREDICTION OF CARTRIDGE DYNAMIC PROPERTIES 

One of the objectives of the work reported is to develop, and 
evaluate, relatively simple predictive methods for practical elastomer 
configurations, such as radially loaded ring cartridges. In this 
Section a number of alternative methods of prediction of radial stiff- 
ness are developed and numerical values are obtained for the particular 
cartridge configurations selected for testing. 

Figure 10 illustrated the geometry of the elastomer ring cartridge. 

The elastomer material separates the inner diameter of one cylinder 
from the outer diameter of a second, smaller cylinder. The loading 
of interest is radial; that is, it tends to move one cylinder radially 
eccentric to the other. 

Three predictive methods have been tested: 

• Beam-Column Method 

• Method of Gobel (ref. 19) 

• Predictions Based Upon Plane Stress Analysis 

Each of these is a ‘static* prediction method; that is, it seeks a 
ratio of element stiffness to storage modulus which is a function of 
the geometry only and applies that ratio in the determination of 
dynamic stiffness and of damping. The approach is not rigorously 
justifiable, but may be given credibility by empirical verification. 

Beam- Column Method 

In this development the cartridge is assumed to be made up of a 
series of noninteracting, elemental column/beams , whose variation of 
thickness due to radius effects is small. Figure 26 illustrates the 
cartridge and the Beam-Column element, whose mean loaded area is 

da = r j6d6 (55) 

where ^ is the length of the cartridge 

r is the mean radius = (r^ + x^)l2. 

Consider the rigid inner member to be displaced radially by a small 

radial distance Y, relative to the rigid outer member. The compression 

deflection of the element is Y cos 0, if 0 is the angle the radial line 

to the element makes with the direction of displacement. The shear 

deflection of the element is Y sin 6. So, considering the element to 

act independently of adjacent elements, the elemental compression and 

shear forces, df , df are 
c s 

df =5 Y cos 0 E 7 ^ ' — r (56) 

c (r^ - r^) 
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Y sin 0 G 


(57) 


df 


r jed9 
(r^ - r^) 


and the total force exerted at the mean radius is 
-,2n 


F = 




Y r jg 


o ^^2 


E cos^ 9 d0 + G sin^ 6 d0 


Y r A tt 


] 


(V2 - rp 


E + G 


(58) 


so that, if we assume E = 3G (Poisson's ratio v = 0.5 implied), then, 

(59) 


K = T7 = 2n ^ Gi BEAM- COLUMN RESULT 

r Y (r„ - r.) 


A slightly different result is obtained if the tapering of beams and 
columns due to radius effects is accounted for. 


K 


4tt Gi ACCOUNTING FOR 


r a (r«/r,) RADIUS EFFECTS 

n z i ■■■■“. 


(60) 


But, as Table IV shows, the reduction in stiffness due to radius effects 
is less than four percent for radius ratios up to two. 


TABLE IV. REDUCTION IN STIFFNESS DUE TO RADIUS EFFECTS 
(Relative to Beam-Column Result) 


r^/ti 

Reduction 

1.1 

0.087. 

1.5 

1.357. 

2.0 

3.827. 


The radius ratio selected for testing is 1.5, for which radius effects 
are only 1.35 percent, and have been neglected. 

Method of Gobel 

The method of Gbbel initially follows a very similar development to the 
Beam-Column method accounting for radius effects, with the intermediate 
result : 


K 

r 


I 


n 


n i 


(G + E) 


(61) 
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However, GObel then claims that, for 3 r = 1, the effective E-modulus 

is given by ^^2 ^ 1 "^ 

Egff = 6.5 G (62) 


and, in addition that, for other length to radius ratios, a correction 
factor, fj^, should be applied, so that 


7.5 r\l Gfj^ 

•'r ■ (rj/rj) 


(63) 


Gobel presents f^ as a graphical function of Payne and 

Scott (ref. 6 ) in their book ENGINEERING DESIGN WITH RUBBER present a 
cubic fit to Gobel 's function: 

f^ = 1 + 0.0097 (64) 

which agrees closely with Gobel *s graphical presentation. The basis 
for or for the statement E = 6.5 G is not made clear in Gobel *s 
book. These are not presented as empirical results, since no test 
data is given for this geometry. Gobel ’s method implies that = 

6.5 G is an asymptotic condition for very short length cartridges. 

This implication is surprising. 

One set of static test data for the cartridge configuration (Adkins 
and Gent, ref. 7) has been identified and this data is compared with 
the prediction of the Beam-Column method and the method of Gobel in 
Figure 27. It shows good agreement at higher values of X/Cr^ ri ) 
between Gobel and test data but at lower ^/(^2 ~ Beam-Column 

method appears to come closer, due to a much sharper drop in stiffness 
of the test element. This comparison suggests that E^ff = 6.5 G is 
not a valid asymptotic condition for short cartridges. 

Predictions Based Upon Plane Stress Analysis 

The value of using a plane stress elastic analysis was also investigated. 
In this case, only a quarter of the cartridge was modeled and appropri- 
ate symmetry and antisymmetry boundary conditions were imposed. The 
quarter cartridge was embedded in a rectangular finite difference mesh 
and a finite displacement of the inner member relative to the outer 
member was imposed. An existing computer program (MTI Program No. 444) 
was used to solve for the internal displacement field within the 
cartridge. Since it is not possible to employ a Poisson’s ratio of 
0.5 in this computer program, a value of 0,4983 was used, together 
with a shear modulus value of 1.795 x 10^ N/m^, This combination of 
Poisson’ s^ratip and shear modulus corresponds to a bulk modulus of 
5.28 X 10^ N/m^, which is within the range of bulk moduli determined 
by Rightmire (ref. 18). 
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REDUCED RADIAL STIFFNESS K/GX 



RATIO OF LEN6H TO RADIAL EXTENT 

^ / (r2“r|) 

I'ij;. 27 Comparison of Prediction Methods for Cylindrical Cartridge 
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The stiffness is calculated by determining the component of shear and 
direct stress in the direction opposing the displacement x and inte- 
grating over the active cartridge surface to obtain the net force. 

The stiffness is, then-, the ratio of net force to displacement. 

Numerical Predictions of Cartridge Stiffness and Damping 

Each of the predictive methods described for the ring cartridge has 
been applied to the geometry selected for testing. The results are 
presented in Table V in nondimensional form; that is, the ratio of 
stiffness to the product of shear modulus and length is presented. 


TABLE V. NONDIMENSIONAL STIFFNESS FOR CARTRIDGES 
SELECTED FOR TESTING (Inner and Outer 
Diameters: 1.905 and 2.858 cm) 


Method 

K/Gi For 1 = 0.476 cm 

K/Gi For JL = 0.953 cm 

Gobel 

58.7 

62.6 

Plane Stress 

35.1 

35.1 

Beam-Column 

31.4 

31.4 


It is seen that only the Gobel method reveals any dependence of this 
nondimensional stiffness on length of the cartridge, and that the level, 
of the Gobel predictions is substantially higher than that of the Plane 
Stress and Beam-Column predictions. 

The assumption necessary to apply the Gobel and Plane Stress predictions 
in the determination of dynamic stiffness and damping is that: 



Gi 


(65) 


where K = Ki +• iK^ is the complex impedance of the cartridge and G' = 
G -I- iG is the complex modulus. That is, the relative distribution of 
stresses is independent of time under sinusoidal loading and displace- 
ment . 


The simple Beam-Column method, without accounting for radius effects, 
since it assumes no interaction between elements, may be developed 
directly using either G or G . In other words, beyond the significant 
assumptions in its development already presented, no further assumptions 
arc necessary in applying it to dynamic stiffness and clamping. 
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Using these assumptions where appropriate, the predictions of Table V 
may be applied to the cartridges tested by the BERM method. The non- 
dimensional stiffness from Table V is simply multiplied by the length 
of the cartridge and the storage and loss modulus values at the fre- 
quency of interest. 

The results of this prediction approach are compared with test results 
later in this report . 
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EXPERIMENTAL DESIGN AND DESCRIPTION OF TESTS 


A forced vibration, resonant mass type of apparatus was used to test 
the dynamic properties of the elastomer samples. The apparatus was 
designed as a base excitation, electromagnetic shaker driven mass- 
spring system which could be brought to resonance for a range of dif- 
ferently sized vibrating masses on top of the elastomer ’spring’. 
Acquisition of data was, however, not limited to the resonance con- 
dition of the system. In fact, data obtained at the resonant frequency 
of each mass-spring combination was just one of several data points 
acquired at several vibration frequencies around resonance, where 
significant amplitude ratios between base and mass vibration existed. 

The following sections present a description of the test rig, the 
design of the elastomer test samples, the instrumentation as it was 
used for the recording of test data, and the range of test parameters 
investigated. 

Description of Elastomer Test Rig 

The basic features of this test rig have largely been developed in the 
course of work conducted under previous phases of NASA Contract NAS3- 
15334. Reference 2 provides a complete description of test rig capabil- 
ities and an analysis of the dynamic characteristics of the test rig 
components. Figure 28 is a schematic showing all the components. 

For the tests reported herein, not all of the available features of the 
test rig were used. In particular, for test frequencies above 100 Hz, 
neither the mass support air spring nor the guide spokes (guide bearings) 
were used. The simplified test rig consisted of a shake table mounted 
base plate to which the elastomer samples were bonded, and a second 
plate bonded to the free side of the elastomer samples. A piston with 
Bellofram seals and various amounts of deadweight were attached to the 
second plate which vibrated freely on the elastomer spring when the base 
plate was subjected to sinusoidal vibrations from the shake table. 

The piston cavity was externally pressurized with air and served to 
preload the elastomer sample in a static sense. An overall view of 
the simplified test rig preload cylinder is given in Figure 29, where 
a small mass is attached to the free end of the preload piston pro- 
truding from the cylinder. Preload was provided on all tests conducted 
with elastomer compression samples and on a few tests of the cartridge 
configuration. (The presence or absence of preload is noted on all 
test results.) All tests with elastomer samples stressed in shear and 
the majority of the cartridge tests were made without static preload 
on the sample, and for these tests, the rig was further simplified by 
removing the piston and seals. Without piston and seal, all test rig 
components and features with interactive or additive dynamic contribu- 
tions to the test system had been eliminated. What remained was a 
dynamically *’clean'* test rig which, of its own, did not contribute any 
identifiable dynamic "noise** or distortion of test results. 
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Fig. 28 Schematic of Elastomer Test Rig Showing All Components 
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The vibration input to the test rig in all tests was obtained from 
a commercially-available electromagnetic shaker system* capable of 
delivering 66,700 N (15,000 lb) force in the sinusoidal mode of 
vibration • 


The ambient temperature within the elastomer test specimen cavity was 
controlled through the combined effect of electric heating cables and 
water-cooled tubes installed in the base plate below the test specimens 
and heated or cooled air blown directly to the test specimen cavity. 

For higher than ambient temperatures, the electric heaters were auto- 
matically controlled by a temperature controller, which obtained the 
necessary feedback from the thermocouple embedded in the metal of the 
base plate. For ambient temperatures down to 2 C, the air supply to 
the test specimen cavity was passed through cooling coils embedded in 
dry ice. 

Details of the temperature control system and of the temperature 
recording system are shown in Figure 30. 

Elastomer Test Samples - Material Selection 

The initial selection of polybutadiene as test material had already 
taken place during an earlier phase of this program, reported in 
reference 2. Restated briefly, the choice of polybutadiene, which is 
classed as Broad Temperature Range (BTR) elastomer, was primarily based 
upon the desire to test a material that was not overly sensitive to 
small variations in ambient temperatures nor to the temperature varia- 
tion induced in the sample through vibration testing. 

The choice of ^lybutadiene among BTR elastomers was based upon the 
manufacturer's recommendation founded on manufacturing considera- 
tions, such as the achievement of consistent material properties from 
batch-to-batch . This material, which carries the manufacturer's desig- 
nation NEX156G, has a nominal hardness of 70 durometers (or Shore A 
Hardness) and has the highest hardness and damping of those available 
from the manufacturer. 

Elastomer Test Samples - Configurations 
Three sample configurations were tested: 

• Shear Type Samples 

• Compression Type Samples 

• Cartridge Type Samples 


Ling Electronics, Model 335A Shaker with PP-35/70 VC Power Amplifier 
SCO-100 Servo Control Center. 

*>V‘^ 

Nichols Engineering, Incorporated, Shelton, Connecticut. 
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Fig. 30 Temperature Control System For Elastomer Tests 






Shear Test Samples . Two shear test samples were used for testing in 
this phase of the contract. One shear test sample was made up from 
eight individual elastomer strips, each measuring 12.7 mm in height, 

48.8 mm in length and 3.2 mm in thickness (0.5 x 1.42 x 0.125 inch) 

(Fig. 31). A second shear sample consisted of four individual elastomer 
strips, each having twice the height (25.4 mm) of the eight-strip 
sample (Fig. 32). 

One individual strip was glued, on its large face to each of the sides 
of a square aluminum block, with the other large face of each strip 
glued to an additional aluminum holding block. A cyanoacrylate glue 
was used (Eastman 910 MHP) . In the eight-strip sample, two such assem- 
blies, each consisting of four strips, were placed on top of each 
other and tested together. In the assembly of the test fixture, the 
four holding blocks were bolted to the shake table, while the elastomer 
supported center block was attached to, and therefore became part of, 
the resonant mass. 

Compression Test Sample . The compression test sample consisted of 
ten individual elastomer cylinders, each measuring 12.7 mm in diameter 
and 6.4 mm in height (0.5 by 0.25 inch). These cylinders were arranged 
in a circular pattern and glued with their faces to two aluminum plates 
(see Fig 33). The bottom plate was anchored to the shake table and the 
top plate connected to the resonant mass. 

Like the shear test sample, the compression test sample had originally 
been made for the investigations reported in reference 2. Details of 
manufacture and test results for two other compression test samples 
made up of three cylinders 3.2 mm in height and of thirty cylinders 
each 12.7 mm in height may also be found in reference 2, 

Cartridge Test Samples . Individual cartridge test specimens were made 
up of an annular elastomer body with rectangular cross-section held 
between the cylindrical bore of a rectangular shell structure and a 
round shaft of smaller diameter. 

The elastomer ring was molded to tVie shaft and bonded to the surround- 
ing shell, which was of split design to facilitate initial assembly and 
later, static preloading prior to testing. 

Figure 34 shows one cartridge specimen before and after assembly of the 
outer shell around the bonded elastomer ring on the shaft. 

Cartridge dimensions were chosen so that four individual cartridges in 
parallel would have a combined spring stiffness that, in combination 
with the minimum mass arrangement in the test rig (top plate only), 
would result in a system resonant frequency of at least 700 Hz. Fre- 
quencies to 1000 Hz could then be achieved by testing to about 1.4 times 
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Fig. 31 Test Assembly of Eight Elastomer Shear Specimens, 
Each 1.27 cm (0.5 In.) High 
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Fig. 32 Test Assembly of Four Elastomer Shear Specimens, 
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Fig. 34 Elastomer Cartridge Specimen Before Bonding of Outer 
Shell (Left) and with Outer Shell in Place (Right) 


the highest resonant frequency and, through mass addition, resonant 
frequencies loT^er than 200 Hz could be achieved. 

A second set of cartridge specimens was designed with identical metal 
dimensions but with only half the elastomer width of the first set. 

For this set the upper limit of the frequency range was reduced in 
proportion to the square root of the stiffness change, or by a factor 
of 1.4. 

The elastomer dimensions (in mm) of the individual specimens for the 
two sets were as follows: 

Outer Inner 

Width Diameter Diameter 

9.5 28.6 19.1 

4.7 28.6 19,1 

A manufacturing drawing of the cartridge elastomer specimens is shown 
in Figure 35. The manufacturing process of the cartridge elements 
evolved from discussions with the manufacturer. Since a split shell 
was desired to allow a positive radial preload to be applied for 
testing, molding of the elastomer to both shaft and outer shell was 
not considered feasible because tensile stresses setup in the elastomer 
due to material shrinkage would not be circumferentially symmetrical. 

In the selected manufacturing process, the elastomer was first molded 
on the shaft and then ground on the outer diameter to the dimension of 
the circular hole in the assembled shell. With a small interference 
tolerance, the outer shell was then clamped around the elastomer and 
the assembly heat cured. Bonding of the elastomer to the shell was 
facilitated by a special coating applied to the bore of the shell prior 
to assembly. (The coating in the shell is clearly visible in Fig. 34.) 

Prior to assembly into a complete test sample, all cartridge elements 
received a uniform preload of 5.3 percent of their uncompressed radial 
heights. Preloading was accomplished by replacing the shims measur- 
ing 1.625 mm between the halves of the outer shell with reduced 
thickness shims only 1.120 mm thick. The assembly of the individual 
elastomer cartridges into the complete test sample is depicted in 
Figures 36 through 38. An assembly of four cartridges was the 
largest number which could effectively be packaged within the exist- 
ing size constraints of the test rig. In Figure 36 the four 
elastomer cartridges are shown next to the top plate to which they 
will be clamped with two clamps each. In the following picture 
(Figure 37) the cartridges have been placed into their locations 
and in Figure 38, the clamps have been placed over the cartridge element 
shafts and locked to the plate with bolts from beneath. At this stage, 
the elevations of the individual cartridge element shells were measured. 
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Fig. 35 Elastomer Test Sample 
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Fig. 37 Elastomer Cartridge Placed in Top Plate (Inverted Position) 
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Fig, 38 


Elastomer Cartridges Assembled 


in Top Plate 
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and shims were prepared for placement between the flat surface of the 
cartridge shell and the base plate (shown standing behind the assembly 
in Fig, 38). With shimming in place, differences in elevations between 
individual cartridge shells were held to less than 0,025 mm. The 
assembly of the second set of cartridge elements (narrow elastomer 
elements) received a very light grinding cut to remove some lack of 
parallelism between the flat shell surfaces and the bottom plate on 
two of the elements. 

Instrumentation 

The uncomplicated measurement requirements of the base excitation 
resonant mass testing method have been expounded before (see reference 
1). Restated briefly, the measurement requirements for the experimen- 
tal investigation of the elastomer dynamic properties are as follows: 

1. Displacement measurement of elastomer support plate 
attached to the vibration table, relative to ground; 

2. Displacement measurement of elastomer support plate 
attached to the resonant mass, relative to ground; 

3. Phase angle measurement between displacement measure- 
ments (1) and (2) above; 

4. Displacement measurement between two elastomer support 
plates (relative amplitude across the elastomer); 

5. Vibration frequency; 

6. Temperature of elastomer support plates and ambient 
air temperature in test specimen cavity; 

7. Vibrational power dissipation across the elastomer 
test sample. 

For convenience, acceleration was measured in (1) and (2) instead of 
displacement, but since the motions were sinusoidal the displacements 
were directly derivable from the accelerations. The very low level of 
motion of the shaker body (from which the shaker table is isolated by 
air springs with a very low resonant frequency) was also measured by 
accelerometers . 

The input acceleration (1) was measured by an accelerometer located on 
the support base for the elastomer test samples. Output acceleration 
on the resonant mass (2) was determined from three accelerometers 
mounted directly to the top plate attached to the elastomer specimens. 

The use of three accelerometers, spaced 90 degrees apart, with their 
output signals displayed on one oscilloscope screen, permitted immediate 
detection of nonaxial motions of the resonant mass. When nonaxial 
motions of the resonant mass occur, they manifest themselves either as 
amplitude or phase angle variations among the three signals, or as a 
combination of both (testing experience revealed that noticeable nonaxial 
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motions occurred only in a few instances at certain test frequencies. 

Test data was generally discarded when nonaxial motions were observed.) 

The displacement measurement (4) between the two elastomer support 
plates (relative amplitude across the elastomer) was accomplished with 
a noncontacting capacitance-type sensor. 

Chromel-alumel type thermocouples were used to measure the temperature 
in the metal directly adjacent to one or two elastomer specimens in 
each assembly at a distance of approximately 1,6 mm from the elastomer 
surfaces. In the shear sample the thermocouples were embedded at the 
midface location. In the compression test samples the thermocouples 
were directly above and below the center point of the flat surfaces of 
the cylindrical specimens. In the cartridge specimens one thermocouple 
was located in the axis of the shaft at the midpoint of the elastomer 
and two additional thermocouples were placed diametrically opposite in 
the metal of the outer shell, approximately 1.6 mm from the outer 
diameter of the elastomer ring and also at the axial midpoint location 
of the elastomer ring. 

Displacement and acceleration data, together with temperatures from 
individual thermocouples were displayed visually in digital form for 
monitoring purposes. As a permanent record, this data was also 
recorded on a 21-column digital printer. 

Amplitude signals, from the four accelerometers and one displacement 
sensor, were sequentially switched by an analog scanner into a two- 
channel tracking filter, which provided a visual readout of vibration 
frequency and of two filtered amplitude signals. The base amplitude 
signal (acceleration) was fed at all times into one of the two channels 
of the tracking filter and from there into a phase meter where it served 
as a reference signal for the measurement of the phase angle of the 
three output acceleration signals and the displacement signal. The 
d-c values proportional to phase angle and amplitude from the phase 
meter and the tracking filter were than converted into digital form in 
digital voltmeters and printed. Each printed line contained phase angle 
and amplitude from one sensor, together with the vibration test fre- 
quency and the temperature indication from one thermocouple. (See 
Figure 39 for a schematic of the data acquisition system.) 

To monitor the power dissipation in the test sample during testing, a 
previously developed power meter was used. This meter utilizes avail- 
able acceleration and displacement signal inputs. The three (unfiltered) 
accelerometer signals from the top plate of the specimen holder (resonant 
mass signals) are averaged, to give a signal proportional to the net 
force on the resonant mass. This signal is therefore almost equal to 


76 


I 



SIGNAL SIGNAL MONITORING SIGNAL PROCESSING 

ACQUISITION AMPLIFICATION AND SWITCHING AND READOUT 
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that portion of the force acting through the elastomer sample which 
is due to dynamic displacement (differing only by dynamic displace- 
raent forces in the air cylinder in those cases where static preload 
is applied to the test sample). The displacement signal, which 
indicates relative motion across the sample, is differentiated to 
give a signal proportional to rate of deformation (relative velocity). 
The force and velocity signals are multiplied together to give the 
instantaneous power. The instantaneous power consists of a d— c 
(average) component and a sinusoidally varying component and, by 
filtering out the sinusoidally varying component, the average power 
dissipation is obtained. 


Testing Procedures 

One of the major features of the current test program was the performance 
of the scheduled tests at near-resonance conditions. Since phase angle 
between base excitation and resonant mass response is an accurate in- 
dicator of the amount of damping in the region of resonance, measure- 
ments should preferably be made in the phase angle range between approx- 
imately 15 and 165 degrees (as discussed in reference 2). This requires 
that the test frequency be approximately 0.7 to 1.5 times the critical 
frequency of the elastomer-resonant mass system. Coverage of the full 
test frequency range of 100 to 1000 Hz with sufficient test points, 
therefore, necessitates changes in the size of the resonant mass. 
Typically, between five and seven different mass arrangements ranging 
from 0,5 to 15 kg were found sufficient to cover the desired test range. 
The smallest mass consisted of the elastomer test sample top plate only, 
with no additional weight attached. Increasingly larger mass combina- 
tions were obtained by attaching weights of up to 15 kg directly to the 
top plate. In the case of the compression sample, which was tested 
with preload applied through the air piston cylinder arrangement, 
weights up to 10 kg were directly attached to the preload piston. 

Once any one of the elastomer test samples had been mounted on the shake 
table and the necessary instrumentation installed and connected, a 
typical test sequence proceeded as follows: 

1. A resonant mass was selected and installed, 

2. In the case of the compression tests, the elastomer 
sample was statically preloaded. 


I 
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3 . 


The elastomer test sample cavity was enclosed and the 
temperature control system given time to adjust ambient 
temperature to the desired value* 

4. With low vibration levels applied to the base of the 
elastomer holding fixture, frequency scans were con- 
ducted until the approximate resonant frequency of the 
system was found, (For selected cases an exact 
determination of the system resonant frequency was 
subsequently made by tracing, on an X-Y plotter, the 
d-c value proportional to the resonant mass accelera- 
tion amplitude as a function of the base excitation 
frequency,) It may be noted here that, for a base- 
excited single-degree of freedom spring-damper-mass 
system, resonance occurs at an angle smaller than 90 
degrees. The deviation from 90 degrees is essentially 
determined by the amount of damping in the system. 

5. While the predetermined vibrational power dissipation 
level in the elastomer test sample was maintained by 
adjustment of the shaker power input level, the vibra- 
tion frequency was reduced until the phase angle de- 
creased to a value between 15 and 25 degrees. 

Provided none of the acceleration and displace- 
ment signals showed signs of abnormalities 
(distortions, or indications of nonaxial motion of the 
resonant mass), a data point was recorded. The recording 
of one data point consisted of one complete printout of 
all amplitudes (up to six individual signals), associated 
phase angles and test frequency, and all temperatures 
from the thermocouples located in the test specimen 
holding plates. 

6. Stepwise increases in vibration frequency were imposed 
and between six and eight data points were recorded 
until the phase angle reached approximately 165 degrees. 

7. Tests, comprising steps 1 through 6, were then repeated 
with each of the remaining masses in turn, each mass 
giving a dynamic system with a different resonant 
frequency. 

8. Shake table power was turned off and, for the compression 
tests the preload on the elastomer sample was increased 
to the next higher value. 

9. Steps 4 through 7 were repeated until data at all desired 
preloads were obtained. 

10. Shaker table power was turned off and, for the compression 
sample, all preload on the elastomer removed. The elastomer 
test sample was then given time (at least one hour) to 
recover its original, uncompressed height (as indicated 
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by displacement probe readings) before tests at either 
higher ambient temperatures or a higher power dissipa- 
tion level were commenced. 


The result of this test sequence is a printed data set defining the 
measured amplitudes for each vibration sensor, and the temperature 
values at each frequency point for a series of mass values. Together 
with data specifying the test conditions, this data set is entered 
directly onto punched cards which are processed by a computer. All 
required editing and averaging functions are performed by the computer, 
then stiffness and damping, , are calculated at each frequency 

as follows: 


2 2 

CD M [g ^ a cos cpj 
2 

a “2a cos cp + 1 


( 66 ) 


K 


2 


2 


a 


CO M g sin cp 
^ 2g cos cp + 1 


(67) 


where a is the transmissibility , that is the ratio of absolute 
values of output to input acceleration; and 
cp is the phase angle between input and output acceleration 
signals . 

The development of these expressions for stiffness and damping is pre- 
sented in Appendix B. 

Test Parameter Ranges 


It is the purpose of this Section to provide an overall guide to the 
test cases reported for this investigation and to indicate the range 
of parameters covered in each case. The method of presentation chosen 
is predominantly that of data plots showing calculated elastomer stiff- 
ness and damping (in N/m) versus test frequency. These plots were 
originally prepared by computer line-printer, and offered a rapid, 
economical means to evaluate the combined effect of such constant 
test parameters as dissipation power level in the specimen and 
ambient temperatures. The purpose of presenting these plots in 
this section is to indicate the extent of parameters, particularly 
in the frequency domain. Lines through individual data points have, 
therefore, not been extrapolated beyond either the lowest or the 
highest frequency data point. (The following section of the report 
presents the test data in expanded format, one frame for each test 
condition, together with theoretical correlations and comparisons. ) 


Tests Performed on Shear Samples 

Tests were performed on two different shear samples. One sample con- 
sisted of eight elastomer strips, each 1.27 cm high and the second 
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sample was made up of four strips each having twice the height 
(2,54 cm) of the eight-strip group. The eight-strip sample was 
tested over a range of ambient temperatures and the four-strip 
sample was tested over a range of power dissipation levels. Test- 
ing of the eight— strip sample was suspended because it was believed 
that chemical decomposition of the elastomer had been initiated 
when the ambient temperature had been allowed to go up to 100-110 C 
during the 93 C tests. Evidence of damage was inconclusive — a 
retest at 32 C (Fig. 40) did not show significantly different be- 
havior from original tests, but subsequent tests at 5 C ambient 
temperature showed elastomer damping behavior that appeared 
anomalous when compared to data obtained at other temperatures. 

The range of test parameters (frequency and ambient temperature) 
covered in tests of the eight-strip sample are shown in Figure 41. 
The ranges of power level and frequency covered in tests of the 
four-strip sample are shown in Figure 42. A listing of shear test 
cases is given in Table VI. 


The power dissipation levels are nominal values, representing the 
average of values calculated, after the fact, from measured frequency, 
damping and amplitude (power ® 0.5 o) Y* ^2). The values for shear 
and compression specimens are generally higher than planned. It was 
intended that each test would be conducted at a constant power level 
but it was identified, towards the end of the series of shear speci- 
men tests, that a constant power meter setting actually resulted in 
a variation with frequency: low frequency dissipation was higher than 

average and high frequency dissipation was lower than average. The 
average value was also higher than intended. This characteristic is 
shown in Figure 43. Following a recalibration of the meter it was 
possible, during the cartridge tests, to maintain more consistent, 
targeted, dissipation levels. 

The highest test frequency is limited by the lightest resonant mass 
without the preload piston. This limit is 0.5 Kg and frequencies for 
the stiff est shear specimen are as high as 1600 Hz. The lowest test 
frequency is governed by the highest mass which can be resonated with- 
out rocking motion and without applying an excessive preload (arbitrar- 
ily limited to one percent). The normal upper mass limit was 15 Kg. 

It may be seen, in Figures 41 and 42 that, with these limitations, the 
test frequency range covered shifts downwards as the stiffness of the 
test specimen falls - either with increasing temperature or with in- 
creasing dissipation. 

Figures 44, 45 and 46 present some interesting comparisons which relate 
mainly to the test method. In Figure 44 test series 4-01-^06 and 
4-19-^24 are compared. The conditions for each set were Identical but 
the two tests were separated by a two week time delay in which the rig 
was disassembled and reassembled. The agreement between the two is 
good and indicates that good repeatability can be achieved with the 
test method. 
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Fig. 40 Stiffness and Damping for Eight-Strip Elastomer Shear 
Specimen at 32 C Ambient Temperature and 0.062 watt/cm 
Power Dissipation 
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TABLE VI. LISTING OF SHEAR SAMPLE TEST CASES 

[Test Conditions Common To All Tests] 
Preload On Sample: None 
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Fig. 45 Stiffness and Damping for Shear Test Specimen at 32 C with 
and Without Lower Air Cylinder of Test Rig 
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Pig. 46 Stiffness and Damping for Shear Specimen at 32 and 66 C 
With Lower Air Cylinder Vent Holes Open and Closed 




In Figure 45 data sets 4-38->-43 and 4-62->67 are compared. The first 
set was obtained without the air cylinder (preload piston) and the 
second set was obtained with the air cylinder in place, and sealed, but 
not pressurized. In other words, apart from the presence of the 
piston, the conditions are identical. The differences between the 
two sets of data are small, but they are undeniably in the direction 
of lower stiffness and damping with the cylinder in place. This 
result is surprising, since analysis of the test rig mechanics in- 
dicates that the air spring stiffness and damping will both tend to 
act in parallel with the test specimen. However, the result was 
fully confirmed by data sets 04-01-t-06 or 04-19-)-24 and 04-31^-35 (not 
shown in graphical presentation). This anamolous trend is not fully 
understood at this time, but is assumed to be the result of flexibil- 
ity in test rig elements (for example, the piston or cylinder head). 

In Figure 46 data sets 4-38-v43 and 4-44-^9 for 32 C are compared and 
data sets 4-50^55 and 4-56->-61 for 66 C are compared. In each com- 
parison the preload piston is present, but between the first and second 
set, the only difference is that vent holes in the cylinder have been 
opened. The effect of this change is clearly negligible, leading to 
the conclusion that anomalous behavior of the air cylinder is not a 
result of pressure fluctuations. 

These anomalies with the preload piston confirm that small uncertainties 
are introduced in tests with the piston and that the most reliable 
elastomer specimen data will be obtained without the preload piston. 

Tests Performed on Compression Samples 

Under Phase III all compression sample tests were performed on an 
assembly of ten elastomer specimens 12.7 mm in diameter and 6.35 mm in 
height (see Fig. 35). Tests were generally conducted at 2--1/2 percent 
preload and then repeated (for each series of masses) at five percent 
preload. At a nominal power dissipation level of 0.062 watt/cm^, the 
compression sample was tested at three ambient temperatures; 2, 66 and 
93 C. At 32 C ambient temperature the sample was tested at Increased 
power dissipation levels of approximately 0.50, 1.1 and 2.2 watt/cm^. 
(Results from tests at 32 C and low nominal power dissipation level 
(0.062 watt /cm ) were available from a previously conducted test series 
in Phase II . ) 

The frequency range covered in compression sample tests extended from 
approximately 150 Hz for high temperature tests to 1000 Hz for low 
temperature tests. Coverage at the high frequency end was generally 
not as high as for the shear sample tests, due to the preload piston 
which increased minimum mass, 

'rest range parameters for all compression sample teats are shown In 
Table VII. 

Graphical representations of the test ranges covered for the compression 
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TABLE VII. LISTING OF COMPRESSION TEST SAMPLE CASES 


Test Case 
Number 

Mass On Elastomer 
(Kg) 

Ambient 

Temper- 

ature 

(C) 


Power 

Dissipation 

Level 

3 

(watt/cm ) 

10-01-05 

0.91, 1.81, 2.78, 3.73, 5.35 

66 

2.5 

0.062 

10-06-09 

0.91, 1.81, 2,78, 3.73, 5.35 

60 

5 

0.062 

10-10-13 

0.91, 1.81, 2.78, 3.73, 5.35 

93 

2.5 

0.062 

10-14-17 

0.91, 1,81, 2.78, 3.73, 5.35 

93 

5 

0.062 

10-18-21 

5.35, 3.72, 1.81, 0.91 

32 

2.5 

0.44 

10-32-36 

0.41, 1.81, 3.73, 5.35 

32 

5 

0.55 

10-25-28 

0.41, 1.81, 3.73, 5.35 

32 

2.5 

1.1- 

10-29-32 

0.41, 1.81, 3.73, 5.35 

32 

5 

1.1 

10-37-39 

0.41, 1.81, 3.73, 5.35 

32 

2.5 

2.2 

10-41-44 

0.41, 1.81, 3.73, 5.35 

2 

2.5 

0.062 

10-45-48 

0.41, 1.81, 3.73, 5.35 

2 

5 

0.062 
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sample are shown in Figures 47 and 48. Both figures are for the 2-1/2 
percent preload test case only; very similar results were obtained 
with five percent preload on the test sample. The consistency of the 
stiffness and damping values obtained for the compression test sample 
appears to be good, following very much the trend lines already observed 
and discussed for the shear test sample. 

Test Performed on Cartridge Samples 

Two cartridge type samples were tested, one with elastomer elements 
9.55 mm wide and a second with elements half as wide (4.75 mm). Both 
samples were first tested at four ambient temperatures: 32, 43, 66 and 
80 C, and then at several power dissipation levels with ambient tempera- 
tures controlled at 32 C. The frequency range covered in these tests 
is indicated, for each ambient temperature and power dissipation level 
test case, in Figures 49 and 50 for the wide cartridge; and in Figures 
52a & 52b for the narrow cartridge. Due to the lower overall stiffness 
of the narrower element sample the test frequency range did shift 
somewhat towards lower frequencies. 

The effect of static preload produced by the air cylinder (2-1/2 and 5 
percent deflection of the height of the elastomer material annulus) was 
investigated for the wide cartridge specimen only. In Figure 51 results 
without the preload piston attached are compared with results obtained 
with the piston attached and pressurized to both 2-1/2 and 5 percent 
preload. The effect of introducing the piston appears to be much stronger 
than the effect of doubling the preload. While no direct comparison, 
showing the effect of the piston alone, is available for the cartridge 
specimen, the indications are, just as for the shear specimens, that 
the introduction of the piston reduces the observed stiffness and damping. 

A listing of all cartridge sample test cases is given in Table VIII. 

Static Tests 

Static stiffness tests were conducted with the elastomer samples 
mounted in the test rig and the displacement probe output used for 
automatic recording of elastomer specimen deflection. Using the same 
environmental control as for dynamic tests, static deflection data 
was obtained at temperatures between 2 and 93 C. Application of load 
was controlled manually by pressurization of the preload cylinder. 

Static deflection tests of elastomer materials differ drastically 
from those associated with most engineering materials. An elastomer 
under load will deflect rapidly to a certain extent as the load is 
applied and then continue deflecting at a much lower rate for an 
appreciable amount of time. The rate of deflection decreases with 
time after each step in load application, but for the presently con- 
sidered material, Increases with increased total load on the sample. 

One example of a typical load deflection relationship, as obtained for 
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Stiffness and Damping for Wide (9.55 am) Cartridge Element 
Specimen at Ambient Temperatures of 32, 43, 66 and 80 C 
and 0.044 watt/cm^ Power Dissipation 
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Fig, 52a Stiffness and Damping for Narrow (4.75 mm) Element 
Cartridge Specimen at 32, 45, 66 and 80 C Ambient 
Temperature and 0.022 watt/cm^ Power Dissipation 





TABLE VIII. LISTING OF CARTRIDGE TEST SAMPLE CASES 



Ambient 


Power 


Temper- 


Dissipation 

Test Case 

ature 

Mass On Elastomer 

Level 

Number 

(C) 

(Kg) 

(watt/cm^) 



Sample Element Width 9.55 mm 


2.06, 2.98, 3.96, 6.51, 11.27, 15.65 


2.26, 3.46, 4.46, 11.77 
2.56, 4.46, 11.77, 16.14 


Sample Element Width 4.75 


4C2-01-«04 

32 

2.06, 3.96, 6.51 

, 15.67 

4C2- 05-08 

43 


4C2-09-12 

66 



4C2-13-16 

80 



4C2- 17-20 

32 



4C2-21-24 

32 



4C2-25-28 

32 



4C2-29-32 

32 



4C2- 33-36 

32 

\ 

f 



0.5 

1.0 

1.8 

0.044 

0.044 

0.044 



*With Piston (2-1/2% Preload). 
**With Piston (5% Preload). 


100 















the elastomer cartridge sample at 66 C, is shown in Figure 53. Beginning 
at time zero with no load on the sample (upper left-hand corner in Figure 
53) the output of the displacement sensor (digital voltmeter reading) is 
high, indicating a large gap between sensor and target. As the air pres- 
sure in the loading cylinder on the test rig is increased, first to 4 
psi and then successively in 4 psi steps to 16 psi, the voltmeter read- 
ing decreases sharply at first and then much more slowly. Approximately 
twelve minutes after application of each load step the plotter pen was 
brought back to the left side of the paper and the pressure in the load- 
ing cylinder was increased. 

From Figure 53 and similar data curves the plots for the shear and cartridge 
samples - Figures 54 and 55 - have been assembled. The load was plotted 
against the last deflection value for each load. Similar loading was ap- 
plied to the compression specimen but, in this case, a mechanical dial 
indicator was used for deflection measurement. Figure 56 shows compression 
specimen deflection data. 

Examination of the load versus deflection curves for the shear sample 

• 54) shows several trends which may also be observed for the com- 
pression sample tests (Fig. 56), but are less pronounced for the car- 
tridge sample (Fig. 55). ^ 

1. The load versus delfection curve is nonlinear at small 
deflections only. 

2. The degree of nonlinearity decreases with decreasing 
elastomer stiffness. 

3. Elastomer stiffness decreases with increasing ambient 
temperature . 

Details of the nonlinear portion of the stiffness curve for the shear 
sample (Pig. 54) were confirmed, with greater precision, through sub- 
stitution of deadweight loading for air cylinder pressurization. 

For the cartridge sample (Fig. 55) the deflection curves do not seem 
to ^straighten out^ within the test range as they do for the other samples. 
The limited range of cartridge sample deflections (about 0.33 mm), may 
be directly responsible, which was the result of displacement sensor 
limitations . 
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Fig. 53 Load vs Deflection Curves for Elastomer Cartridge 
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FIr. 56 Static Deflection va Force for Compression Element Specimen 
at Various Temperatures 
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TEST RESULTS 


The tests described in the previous sections generated discrete values 
for stiffness and damping of shear, cylindrical and ring cartridge 
elastomer elements at different values of frequency, temperature and 
dissipation level. The purpose of this section is to present this 
data in its entirety, as well as to provide an appropriate summary of 
it. The data presented forms the basis for further interpretation. 

The basic medium of presentation is the plot of stiffness and damping 
versus frequency. One frame is presented for each set of conditions 
of: 

• Geometry 

• Temperature 

• Dissipation Level 

• Preload (where appropriate) 

On each frame discrete stiffness and damping values are plotted, as 
well as a line representing the least squares best fit power law re- 
lationship between stiffness, damping and frequency, in the form: 

CD (68) 

®2 

= A 2 CD ^ (69) 

The shear specimen power law expressions for different temperatures 
are used to develop a catalog of material properties for this material. 
The material properties catalog is subsequently used as part of the 
basis for predicting stiffness and damping of the cylindrical compres- 
sion and cartridge specimens. 

A brief description of the various sets of test results follows. 

Shear Specimen Tests at Different Temperatures 

These tests were performed with the eight-element shear specimens, 
each element having a dimension of K27 cm in the direction of loading. 
The tests were performed with the simplest configuration possible - a 
mass symmetrically supported by the shear elements which are mounted 
on the table. The preload piston which has been used in previous 
tests was not used, eliminating any effects of tare damping. 

The data is presented in Figures 57 through 60 for temperatures of 32, 
50, 66 and 93 C. The consistent effect of temperature is to reduce 
stiffness and damping. The data shows reasonable consistency for the 
three higher temperatures, but at 32 C some scatter is encountered at 
the highest frequency. The power law relationship provides, in each 
case, a reasonable fit to the data and, certainly, no more sophisticated 
representation would be justified. 
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Fig. 58 Shear Specimen at Zero Preload, 50 C (Eight^Elements as 
Shown). Nominal Dissipation 0.062 watts/cm 
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Fig. 59 Shear Specimen at Zerp Preload, 66 C (Eight Elements as 
Shown). Nominal MsslpatAnn 0.062 watts/cm^ 
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Fig, 60 Shear Specimen at Zero Preload, 93 C (Eight^Elements as 
Shown), Nominal Dissipation 0.062 watts/cm 
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Material Properties Catalog 
Applying to the shear specimen 


G' 


G" 

where h is the strained dimension (0.3175 cm), 
nA is the total stressed area (49.55 cm^), 

A power law relationship between storage and loss moduli (g', (f') and 
frequency has been obtained, as shown graphically in Figures 61 and 62. 
The coefficient values are shown in Table IX. 
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TABLE IX. MATERIAL PROPERTIES CATALOG FOR 
POLYBUTADIENE, NICHOLS NEX 156G 


Temperature 

Storage Modulus 

Loss Modulus 

32 C 
43 C 
66 C 
93 C 

r' - -i c.Q£. ir.6 0.2037 

G = 3.686 X 10 0) 

G' = 1.904 X 10^ 

G' = 1.097 X 10^ 

G" = 1.197 X 10^ 

— Q 0 0 Q 1 ^ • 1^77 

G = 0.333 X 10 CD 

-5 07 1 in^ —0.06372 

G = 3 ,971 X 10 CD 

7/00 1 *“0 .03011 

G - 2 .423 X 10 CD 

1 +0.2246 

G = 1 .064 X 10 CD 


Together with the geometry factors developed in preceding sections, 
this material properties data has been applied in the prediction of 
stiffness and damping of the cylindrical and cartridge specimens. For 
simplicity of application, it is convenient to have available the 
values of G , G at 100 and 1000 Hz, and these values, based on Table 
IX are presented in Table X. 


test data the interpretive formula: 

(70) 

(71) 


h 

nA 


K 2 h 


nA 


TABLE X. VALUES OF STORAGE AND LOSS MODULI FOR 
POLYBUTADIENE, NICHOLS NEX 156G AT 
100 AND 1000 HERTZ 



G' 100 Hz 
N/m2 

G' 1000 Hz 
N/m^ 

g" 100 Hz 
N/m^ 

G" 1000 Hz 
N/m^ 

32 C 
43 C 
66 C 
93 C 

1.369 X 10^ 
1.035 
0.7755 
0.5803 \l' 

2.189 X 10^ 
1.894 
1.5606 
1.0201 i 

3.659 X 10^ 

2.633 

1.996 

1.2297 vl 

2.728 X 10^ 

2,274 

1.862 

1.2949 J, 
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Fig. 61 Influence of Temperature on Storage Modulus (Shear Specimen 
Data, Eight Elements) 
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Compression Specimen Tests at Different Temperatures 

This series of tests was performed with the ten-element compression 
specimens, each element 1.27 cm in diameter and 0.635 cm high. In 
this case, the preload piston was a necessary feature of the tests to 
satisfy the requirement that the compression elements not be allowed 
to experience tension. Geometrical preload values of 2.5 and 5 per- 
cent were applied. Temperature values tested were 66, 93 and 5 C. 

These test results are presented in Figures 63 through 68. Again, 
reasonable consistency between the test data and the fitted line was 
obtained. The consistent trend, as for the shear specimens, was that 
stiffness and damping fall with increasing temperature. 

Included in the figures for 66 and 93 C are predictions of stiffness 
as a function of frequency using the empirical geometry coefficients 
developed on the basis of previous (Phase II) test data and using 
the material property values obtained from the present shear specimen 
tests, as presented in Table IX. In addition, predictions using the 
static shape coefficient value (p ^ = 0 =2) are presented on this 

figures . 

Very similar predictions are obtained with either approach. The pre- 
dictions of stiffness are within ten percent of measurement at 66 C and 
within twenty percent at 93 C. Predictions of damping are within 25 per- 
cent at 66 C and within twenty percent at 93 C. No predictions are made 
at 5 C since material property data was not generated at this temperature 

While the stiffness predictions show similar trends to the measurements, 
in terms of rate of increase with frequency, the damping predictions in- 
dicate a slight decrease with frequency and the measurements reveal a 
more pronounced increase with frequency. 

Compression Specimen Tests at Different Dissipation Rates 

These tests were performed with the ten-element compression specimens 
at 32 C. Preloads of both 2.5 and 5 percent were investigated by using 
the preload piston. 

The variation of stiffness and damping with frequency at different 
nominal dissipation levels is presented in Figures 69 through 73. The 
corresponding variation of amplitude is shown in Figures 74 and 75 and, 
using these two sets of data, the details of each test point in terms 
of amplitude, frequency and temperature may be recreated. 

The general effect of increasing dissipation is to decrease both stiff- 
ness and damping, at all frequencies, for these compression specimens. 

As previously discussed, it is hypothesized that a major reason for 
this degradation is the increase in the material temperature. The 
comparison between this measured data and predictions of the simple 
thermoviscoelastic model will be presented in the next Section of the 
report. 
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STIFFNESS AND DAMPING, N/m 



Fig. 63 Stiffness and Damping vs Frequency. Compression Specimen 
(Ten Elements with Dimensions as Shown). 2-1/2 Percent 
Preload, 66 C, Nominal Dissipation 0.062 watts/cm^. 
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FREQUENCY, Hj 


Fig. 64 Stiffness and Damping vs Frequency. Compression Specimen. 

Five Percent Preload « 66 C. (Ten Elements with Dimensions 
as Shown). Nominal Dissipation 0.062 watts/cm^. 2-1/2 
Percent Preload. 
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Fig. 65 Stiffness and Damping vs Frequency. Compression Specimen. 

2-1/2 Percent Preload, 93 C. Nominal Dissipation 0.062 
watts/cm^ 
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STIFFNESS AND DAMPING, N/m 



FREQUENCY, Hz 


Stiffness and Damping vs Frequency. Compression Specimen 
Five Percent Preload, 93 C. Nominal Dissipation 0.062 
watts/cm"^ 




STIFFNESS AND DAMPING, N/m 



Fig, 67 Stiffness and Damping vs Frequency. Compression Specimen. 
5 C. 2*«l/2 Percent Preload; Nominal Dissipation 0.062 
watt/ cm^ . 
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STIFFNESS AND DAMPING, N/m 



Fig. 68 Stlffnes. and Damping va. Frequency. Compression Specimen. 

5 C, 5 Percent Preload; Nominal Dissipation 0.062 watt/cm^. 
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FREQUENCY, HZ 


Fig. 69 Stiffness and Damping vs Frequency. Compression Specimen. 
32 C. Nominal Dissipation 0.44 watt/cm^ 
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FREQUENCY, HZ 

Fig. 70 Stiffness and Damping vs Frequency. Compression Specimen. 

2-1/2 Percent Preload, 32 C. Nominal Dissipation 1.1 watts/cm 
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EEPRODUCmiLITY OF THE 
ORIGINAL PAGE IS POOR 



600 800 1000 


Damping vs Frequency 
rcent Preload, 32 C. 






STIFFNESS AND DAMPING, N/m 



Fig. 72 Stiffness and Damping vs Frequency. Compression Specimen. 

5 Percent Preloadi 32 C. Nominal Dissipation 0.55 watts/cm^ 


STIFFNESS AND DAMPING. N/m 



FREQUENCY, HZ 


Fig. 73 Stiffness and Damping vs Frequency. Compression Specimen. 

5 Percent Preload, 32 C. Nominal Dissipation 1.1 watts/cm^ 



AMPLITUDE, METERS (PEAK-TO-PEAK) 



fig. 74 Asqplltnde as a Function of Frequency. Compression Specimen. 

Pfircent Preload » 32 C (Ten Elements with Dimensions 
as Shovm) 
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AMPLITUDE , METERS ( PEAK -TO- PEAK) 



Fig. 75 Amplitude as a Function of Frequency. Compression Specimen. 

5 Percent Preload, 32 C (Ten Elements with Dimensions as Shown). 
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Shear Specimen Tests For Different Dissipation Rates 


These tests were performed with four-element shear specimens. Each 
element is 2.54 cm long in the direction of loading, the strained di- 
mension is 0.3175 cm thick, and the width of the elements is 4.877 cm. 

As with the eight-element shear specimen, it is possible to perform 
these tests without the uncertainties introduced by the presence of the 
preload piston. All the high dissipation tests were performed at a 
temperature of 32 C. 

As a reference, a low dissipation level of 0.035 watt/cm has been 
tested, followed by nominal dissipation values 5, 10, 20, 50 and 125 
times the reference. The data, in terras of stiffness and damping, is 
presented in Figures 76 through 81 and the amplitudes as a function of 
frequency in Figure 82. 

Reasonable consistency between the test data and the fitted line is ob- 
tained, with scatter of the order of ten percent of the stiffness and 
damping level. Scatter increases at the highest dissipation levels. 

It is generally a trend that both stiffness and damping fall with in- 
creasing dissipation level. 

The slope of the lines also tend to increase with increasing dissipation 
levels - indicating a more pronounced effect of high dissipation at low 
frequencies. This could result from the power meter characteristic dis- 
cussed in the previous Section. 

Cartridge Tests at Different Temperatures 

These tests were preformed on two different cartridge test specimens. 

Each consisted of four elements mounted in the same fixtures and holders. 
The inner and outer diameters of the elements were the same for the two 
specimens (1.90 cm and 2.86 cm); the difference was the length - 0.47 cm 
in one case and 0.95 cm in the other. Data for the shorter elements are 
presented in Figures 83 through 86 and for the larger elements in Figures 
87 through 90, The temperature range covered is 32 to 80 C. The energy 
dissipation rate is 0.022 watts/cm^. 

The data is generally consistent with the fitted power law line, except 
that at low frequencies, for the shorter length elements, the damping 
values fall with frequency before starting to increase with frequency 
above 150 to 200 Hz. The influence of increasing temperature is to de- 
crease the stiffness and damping. In general, the influence of frequency 
is to increase both stiffness and damping at a similar rate. 

Superimposed on these test data plots are predictions using the various 
methods of prediction for the cartridge geometry previously discussed . 

In all cases, the Gbbel predictions and Beam-Column predictions are in- 
cluded. In most cases the plane stress line is also Included. 
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Fig. 76 Stiffness and Damping vs Frequency. Shear Sneclmens 

Zero Preload, 32 C. Nominal Dissipation 0.035 watts/cm^ 
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Fig. 77 Stiffness and Damping vs Frequency. Shear Specimen. 

Zero Preload, 32 C, Nominal Dissipation 0.176 watts/cm' 
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Fig. 78 Stiffness and Damping vs Frequency. Shear Specimens. 

Zero Preload, 32 C. Nominal Dissipation 0.352 watts/cm' 
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Stiffness and Damping vs Frequency. Shear Specimens. 
Zero Preload, 32 C. Nominal Dissipation 0.704 watts/c 
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Fig. 80 Stiffness and Damping vs Frequency. Shear Specimens. 

Zero Preload, 32 C. Nominal Dissipation 1.76 watts/cm' 
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Tig. 82 Amplitude as a Punction of Frequency. Shear Specimen 
Zero Preload, 32 C 





Fig. 83 Stiffness and Damping for Cartridge Specimen. 32 
Nominal Dissipation 0.022 watts/cm*^ 
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Fig. 84 Stiffness and Damping for Cartridge Specimen. 
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Fig. 85 Stiffness and Damping for Cartridge Specimen. 66 C. 
Nominal Dissipation 0.022 watts/cm^ 
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STIFFNESS AND DAMPING. N/m 


PREDICTION 



1 / n Ar\ 1057 

K I = 8.40 X 10 u 


BEAM COLUMN 
PREDICTION 


QOBEL PREDICTION 


PLANE STRESS PREDICTION 


^ PLANE STRESS — 

PREDICTION I 

BEAM COLUMN PREDICTION 


EUSTOMER 


0.9525 cm 


m 


FREQUENCY , HZ 


Fig. 87 Stiffness and Damping for Cartridge Specimen 
Nominal Dissipation 0.022 watts/cm^ 
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Fig. 88 Dynamic Stiffness and Damping for Cartridge Specimen. 44 C. 

Nominal Dissipation 0.022 watts/cm^ (Four Elements as Shown). 
(D^ensions as Shown) . 
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Fig. 90 Dynamic Stiffness and Damping for Cartridge Specimen. 80 C. 

Nominal Dissipation 0.022 watts/cm^ (Four Elements as Shown)- 



soni5 ©xcGptxonSj fox dsniping of tliG sfiortGx spccinienj thG fosti 
data lies within the band defined by the highest and lowest predic- 
tion lines. The shorter specimen data generally approaches the 
Beam-Column predictions most closely, while the longer specimen data 
more generally lies midway between the Beam-Column and Gobel predic- 
t ion lines . 

The stiffness data shows a similar slope to the predictions, whereas 
the damping test data shows a different trend from the predictions. 

The material properties for the elastomers obtained from shear speci- 
men tests indicate a slight decrease in loss modulus with frequency. 
Since the cartridge damping predictions are related to the loss modu- 
lus via a frequency independent geometry factor they also reflect this 
decrease with frequency. The cartridge test data, however, shows a 
consistent increase with frequency. 

Cartridge Tests at Different Dissipation Levels 

These tests are performed with the same two cartridge test specimens 
as for the tests at different temperatures. The temperature for these 
tests was 32 C. Data for the shorter specimen are presented in Figures 
91 through 96 and for the larger specimen in Figures 97 through 101. 

In Figures 96 and 101, corresponding to each plot of stiffness and damp- 
ing versus frequency, is a line showing amplitude as a function of 
frequency. With one exception, there is good distinction between the 
lines on Figures 96 and 101. However, in Figure 96, above 210 Hz, 
the points identified as 1.32 watts/cm^ nominal dissipation lie very 
close to the line for 0.88 watts/cm^. Apparently the required test 
conditions for 1.32 watts/cm^ were not maintained above 210 Hz. Com- 
parison of Figures 94 and 95 confirm that, below 210 Hz, there is sig- 
difference in the stiffness and damping data, whereas, above 
210 Hz, the two figures present very similar results. Figure 95 is 
intended, for completeness since, toghether with the data of Figure 
96 it fully defines the test conditions for individual points. 

The effects of dissipation level are, generally, to decrease both stiff- 
ness and damping although this trend is less pronounced at high fre- 
quency than at low frequency. In other words, the lines of stiffness 
and damping versus frequency becomes steeper at high dissipation levels. 

Specimen Catalog 

The power law curve fit coefficients for all the specimens whose stiff- 
ness and damping are plotted in Figures 57 through 101 are summarized 
as a function of temperature, dissipation and geometry in Tables XI, 

XII and XIII. These tables form the Specimen Catalog for shear, com- 
pression and cartridge geometries. * 

Transient Behavior 

In high dissipation tests with the compression specimens, a level of 
dissipation was reached at which it was no longer possible to maintain 
a subresonant condition. When the input acceleration to the shake table 
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Fig. 91 Stiffness and Damping for Cartridge Specimen. 
Nominal Dissipation 0.11 watts/cm^ 
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Fig. 92 Stiffness and Damping for Cartridge Specimen, 32 
Nominal Dissipation 0,22 watts/cm^. 
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Fig. 94 Stiffness and Damping for Cartridge Specimen. 32 
Nominal Dissipation 0.88 watts/cm^ 







AMPLITUDE, METERS 



Fig. 96 Amplitude vs Frequency. Narrow Cartridge Specimen. High 

Dissipation Tests, 32 C (Four Elements with Detail as Shown) 


150 


I 



6x10 



u>/N ‘SNIdNVa QNV SS3NddliS 

! 


151 


Fig. 97 Dynamic Stiffness and Damping for Cartridge Specimen. 32 C 
Nominal Dissipation 0.16 watts/cm^ (Four Elements as Shown) 
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>0 200 300 400 600 8001000 2000 4000 6000 10,000 

FREQUENCY, HZ 

Fig. 99 Dynamic Stiffness and Damping for Cartridge Specimen. 32 C 
Nominal Dissipation 1 .00 watts/cm^ (Four Elements as Shown).' 
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Fig. 100 Dynamic Stiffness and Damping for Cartridge Specimen, 32 
Nominal Dissipation 1.8 watts/cm^ (Four Elements as Shovn) 



AMPLITUDE, METERS 



FREQUENCY, HZ 

AmpLiHidi' vs Fi'i-’quency . Wide C;irLridge Specinu'n. High 
Dissipation Tests. 32 C (Four Elements with Detail as Shown) 
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TABLE XI. SPECIMEN CATALOG 
(SHEAR SPECIMEN) 


















TABLE XII. SPECIMEN CATALOG 

(COMPRESSION ELEMENTS) 
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reached a critical level, a slow but consistent change in system 
conditions would set in. With no further adjustments to frequency 
of input acceleration, the output acceleration would steadily increase 
with time, reach a peak, then reduce and settle at a steady value. 
During this change iri acceleration the phase angle would increase 
from as low as 30^ and stabilize at a value above 90^. This change 
in phase angle was a clear indication that the operating condition of 
the system had changed from subresonant to superresonant . This tran- 
sient phenomenon was also accompanied by a significant increase by as 
much as 10 C in the temperature measured near the specimen. 

With a higher input acceleration, the same phenonmenon would occur 
but its time scale would be reduced. When just above the critical 
acceleration level, the peak amplitude would occur after several min- 
utes. However, levels of acceleration could be reached where the peak 
amplitude would occur within 15 seconds. 

The phenomenon was not observed in the shear specimen tests or in the 
cartridge specimen tests. 

Some typical traces of output acceleration as a function of time are 
shown in Figure 102, The frequency in all cases is 433 Hz and the 
supported mass is 0.908 kg. The lowest line on the figure illustrates 
the result of an input acceleration below the critical level at which 
the trnasient phenomenon sets in. In the center line on this figure 
the input acceleration is just above the critical level. Over several 
minutes the output acceleration increases at a slow steady rate and 
then increases sharply until a peak is reached, followed by a gradual 
drop to an asymptotic value slightly lower than the peak value. In 
the upper line of this figure, the rate of increase in output acceler- 
ation is high from outset. The peak is reached within 15 seconds but 
the gradual drop in amplitude to an aymptotic value still takes over 
a minute. 

Phase angles at the beginning and end of each trace are shown as cir- 
cled values. The increase in phase angle from below to above 90® for 
the two lines with a peak is apparent. 

The fact that. this is a transient phenomenon, the fact that a chnage 
in temperature is involved, and the fact that the time scale of the 
phenomenon is controlled by the input acceleration and hence by the 
dissipation rate, strongly suggests that the cuase is a time-varying 
change in the thermal state of the elastomer induced by internal 
heating . 
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MASS ACCELERATION, g/l.4l 



TIME, MINUTES 

Fig. 102 Thermal Transient. Compression Specimen. Supported Mass, 
0.908 Kg, Frequency 433 Hz 
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DISCUSSION 


High Dlaeipatlon Results 

As shown by the test data of the preceding section, the influence of 
a high dissipation rate is generally to reduce the effective stiffness 
and damping of an elastomer specimen. An obvious hypothesis is that 
internal heating reduces local material properties, and, on the, basis 
of this hypothesis, the thermoviscoelastic model described earlier was 
developed. The hypothesis can now be tested by comparing the predic- 
tions of this model with the test data. 

A parameter required by the model is thermal conductivity, and no value 
of conductivity is available for the material tested. However, since 
all other parameters of the model can be assigned values, once material 
properties and geometry factors are available, the approach taken has 
been to exercise the model for a range of reasonable values of con- 
ductivity, and to determine the value giving the best agreement. The 
initial object of this comparison is the compression test data. In 
Figures 103, 104 and 105 stiffness and damping are plotted against 
amplitude for frequencies of 175, 380 and 600 Hz, Superimposed on the 
compression specimen test data are predictions of the model for conduc- 
tivity values of 20 x 10"^ and 6.7 x 10"^ watts/(cm C). It may be 
seen that the predictions for a conductivity of 6.7 x 10"^ match the 
test data with reasonable success, both in terms of trends and in 
absolute values of stiffness and damping. 

The value of 6.7 x 10"^ for conductivity is not unreasonable (rubber 
conductivity is 1.49 x 10"^) but there is no evidence, apart from the 
empirical comparison that this value is correct for polybutadiene nor, 
therefore, is this comparison absolute evidence of the validity of the 
model . 

The observed transient change in amplitude and phase angle under cer- 
tain high dissipation conditions, as presented in the Results Section, 
is qualitatively predictable by the one-dimensional thermoviscoelastic 
analysis described earlier. As an exercise for the model, an attempt 
was made to reproduce the conditions of Figure 102, and to predict 
comparable transient results. In addition to conductivity no value is 
available for specific heat for this material, although a density value 
was obtained from the manufacturer (0.043 Ib/in.^). The conductivity 
and specific heat values used were 6.7 x 10“^ watts/(cm C) and 0.0622 
nm/kg C. The frequency was 433 Hz and the supported mass was 0.908 kg. 

Figure 106 presents some of the results of this prediction for two 
levels of input acceleration (17.25 and 24 g), which are similar to 
those applied to produce the transient phenomenon in Fig. 102. 

A thermal transient phenomenon with a transition through resonance 
is predicted for both input accelerations, but the resonant peaks are 
much sharper than those found in a similar time scale in the corre- 
sponding test conditions shown in Figure 102. The predicted internal 
peak temperatures are surprisingly high, although no measurements of 
internal temperature are available for comparison. 


161 




ui/N'9NldNVa QNV SSBNddliS 


162 


Fig, 103 Stiffness and Damping vs Amplitude for Compression Specimen, 
175 Hz. Comparison of Prediction and Measurements for Two 
Different Values of Thermal Conductivity 
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Fig. 104 Stiffness and Damping vs Amplitude for Compression Specimen, 
380 Hz, 32 C. Comparison of Prediction and Measurement for 
Two Different Values of Thermal Conductivity 
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Fig. 105 Stiffness and Damping vs Amplitude for Compression Specimen, 
600 Hz, 32 C. Comparison of Prediction and Measurement for* 
Two Different Values of Thermal Conductivity. 
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Fig, 106 Prediction of Thermal Transient for ' Comparison With 
Figure 102, 433 H», 0.908 Kg, Supported Mass. Com- 
pression Specimen. Conductivity 6.7 x 10~^ watts/cm C 
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The inability to completely recreate the transient conditions is not 
surprising, since very close simulation of the mechanical properties 
of the elastomer specimen and its dependence on temperature and dis- 
sipation would be required. However, the ability qualitatively to 
predict the phenomena on the basis of a thermoviscoelastic model with 
similar levels of input acceleration is evidence of the importance of 
thermal effects in determining the high dissipation behavior of the 
compression specimens. 


A further check on the model was sought by comparison of its predic- 
tions with test data for the shear specimen. The same value of con- 
ductivity most successful for the compression specimens has been 
applied (K = 6.7 x 10"'^ watts /cm C) . Figures 107 through 110 show 
typical results which exhibit significant differences, both qualitative 
and quantitative, between tests and predictions. The predictions give 
much higher stiffness than the test data. It turns out that the area 
for heat flow is so large and the length of the heat flow path so 
small in the shear specimens, that temperature levels are much lower 
than in the compression specimens. Thus, the predictions based on 
internal heating and a corresponding reduction of material property 
values show only minor deviations from the isothermal values. 


It becomes clear that some other mechanism is influential in reducing 
stiffness and damping of the shear specimens. Reference to available 
literature on the subject shows that another mechanism has, indeed, 
been recognized though not, as yet, analyzed or understood. This 
mechanism shows itself as a softening of shear elements at moderate 
to high strains. A transition region in amplitude has been identified 
as lying between strains of 0.001 and 0.005. Below this strain level 
the stiffness and damping are independent of strain. Above this level 
of strain a fall-off in stiffness occurs. It is reported by Warnaka 
(ref. 4) that consistent behavior is obtained for different materials. 


This softening effect has not been previously reported at these fre- 
quencies for compression loading and there is apparently some dif- 
ference in the behavior of compression and shear specimens with respect 
to high dissipation. Thermal effects clearly account for more of the 
stiffness degradation in compression elements than in shear elements. 


In an attempt to clarify the relative importance of thermal effects 
and the ' strain softening ' effect, stiffness and damping data for the 
compression and shear specimens has been plotted, separately, against 
strain and against dissipation rate, as independent variables. The 
results are shown in Figures 111 and 112 for shear, 113 and 114 for 
compression. The compression specimen plots include data from the 
Phase II series of tests. Some interesting observations may be made: 


1. The shear specimen stiffness and damping data shows 
reasonably close correlation with strain. 
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Comparison of Predictions and Measurements of Stiffness and 
Damping Under High Dissipation Conditions. Shear Specimen. 
(Four Elements) 




168 


Jig. 108 Compaxison of Predictions and Measurements of Stiffness and 
Damping Dnder High Dissipation Conditions. Shear Specimen. 
(Four Elements) . 
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Fig. 109 Comparison of Predictions and Measurements of Stiffness and 
Damping Under High Dissipation Conditions. Shear Specimen. 
(Four Elements) 



Fig. 110 Comparison of Predictions and Measurements of Stiffness and 
Damping Under High Dissipation Conditions. Shear Specimen. 
(Four Elements) 
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STRAIt^, m/m 

Fig. Ill Stiffness and Damping as a Function of Strain. Shear 
Specimen. 32 C (Four Elements) 
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DISSIPATION, WATTS /CM^ 

Fig. 112 Stiffness and Damping as a Function of Dissipation. Shear 
Specimen. 32 C (Four Elements) 
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STRAIN (DOUBLE AMPLITUDE), m/m 

Stiffness and Damping as a Function of Strain. Compression 
Specimen. 32 C. (Ten Elements) 
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STIFFNESS AND DAMPING, N/m 






2. The sheer specimen data shows poor correlation with’ 
dissipation rate. 

3. The compression specimen data shows reasonably close 
correlation either with strain or with dissipation 
rates . 

4. The data is generally at a high enough strain rate 
that strain independence is not apparent, (although 
more extensive data would be required, at low 
strains, to determine where strain independence does 
occur) . 

5. Damping values exhibit a peak, in the case of the 
shear and compression specimens, for strains in the 
region of one percent. 

In summary, these results justify the following statements; 

Stiffness and damping of elastomer elements in the frequency range 
100 to 1000 Hz are influenced by two aspects of vibration amplitude: 
firstly strain, which has a ^softening effect’ above a certain transi- 
tion value; secondly, self-heating due to thermal dissipation, which 
causes reduction in values of local material properties. The rela- 
tive importance of these two effects is governed by the geometry, 
type of loading, and, probably, by frequency and material. The strain 
softening effect appears to dominate for the shear specimens of poly- 
buta diene . 

Further information will be required before the relative contribu- 
tions of thermal and strain softening can be established for the 
compression specimens. The conductivity must be separately deter- 
mined before the one-dimensional thermoviscoelastic model can be 
tested further. The low strain region must be more extensively 
explored to verify the existence of a strain ’transition* region, to 
determine what conditions are required to establish strain indepen- 
dence, and to establish the influence of geometry on the strain 
dependence. Temperature measurements within the material under high 
strains would provide additional insight. 

The present indications are that high amplitude stiffness and damping 
follow closely a unique relationship with strain. If this observa- 
tion can be supported with further test data it offers an interesting 
and simple approach to accounting for high amplitude effects in 
design. 

Damping Performance 

The effect of strain on damping performance is further illustrated in 
Figure 115. Here the loss coefficient, T], has been plotted as a func- 
tion of strain. The loss coefficient is a good measure of the 
elastomer’s ability to damp vibrations of a mechanical system and, 
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Fig. 115 Loss Coefficient as a Function of Strain. Elastomer Specimen 



for a one-degree-of-freedom system such as is used for the Base Ex- 
citation Resonant Mass tests, the damping ratio {Q = C/C ) is equal 
to one-half of the loss coefficient. Consider Figure 115: for 

strains below 0.02 there is a universal increase in loss coefficient 
for all specimens; at a strain of 0,02 all specimens show a loss co- 
efficient in the range 0,21 to 0.26 (10,5 to 13 percent of critical 
damping). Above a strain of 0.02, however, there is some divergence 
in loss coefficient behavior between the different specimens. The 
shear specimen reveals a steady increase in loss coefficient to a 
value of 0.3 at a strain of 0.1. The cartridge specimen loss coef- 
ficient remains constant and the compression specimen loss coefficient 
falls above a strain value of 0.02 

Based on these results, the generally best operating condition is with 
a double amplitude strain of 0.02. Of course, the question of surviv- 
^bility and life must be assessed further before any overall conclu- 
sion about operating conditions are drawn. For shear elements the 
temperature developed due to self-heating are unlikely to be a problem, 
but fatigue life for significant strain amplitude at frequencies in 
the range 100 to 1000 Hz must be determined to be acceptable. For 
compression elements and probably for cartridge elements, the temper- 
atures developed are likely to be a problem in addition to fatigue 
life. 

It is of particular interest that the two different cartridge geom- 
etries show such similarity of loss coefficient dependence on strain. 

It is also of importance to note that the cartridge geometry indicates 
the least overall sensitivity of loss coefficient to strain, and 
achieves a loss coefficient value of 0.2 or greater for all strains 
above 0.003. 

High Temperature Results 

As expected the test results show a significant influence of tempera- 
ture on stiffness and damping and on the material properties inferred 
from the shear specimen results. The method of reduced variables (the 
WLF law) was discussed in a previous section of this report as a means 
of unifying the effects of temperature and frequency. 

The availability of test results for material properties over a range 
of temperatures now allows the method of reduced variables to be tested 
for this material. 

As identified the method of reduced variables requires knowledge of 
the characteristic temperature. No information is available as to the 
characteristic temperature or transition temperature for the material 
tested and it has been necessary to seek to infer the characteristic 
temperature from the test data. To do this the method has been applied 
for a range of characteristic temperatures. Reduced storage and loss 
moduli, referred to 32 C, have been calculated for each test data point 
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for each temperature tested. The result is a set of discrete, reduced 
storage and loss modulus values, covering a wide range of reduced fre- 
quency values (from 11 to 1600 Hz), referred to 32 C. Power law re- 
lationships between reduced moduli and reduced frequency have been 
fitted to the data and, as a measure of the consistency of the reduced 
data, the rms deviation between test data and the power law lines have 
been calculated. The characteristic temperature has then been sought 
as that value which minimizes the rms deviation. 

In Figures 116 and 117 rms deviation, for reduced storage and loss 
moduli, respectively, are plotted against characteristic temperature. 

The deviation for storage modulus shows a minimum at a characteristic 
temperature of -50 C, The deviation for loss modulus however shows 
no sign of possessing a minimum. There is a continual decrease in 
deviation as the characteristic temperature is increased from -100 to 
+32 C. 

Accepting the storage modulus minimum as indicating the transition 
temperature and the loss modulus data as anomalous, reduced moduli 
have been plotted against reduced frequency for a characteristic tem- 
perature of -50 C in Figures 118 and 119. These plots reflect the 
minimization results. The reduced storage modulus follows, with 
reasonable consistency a single straight line on a log- log plot; the 
reduced loss modulus data remains a series of individual curves 
corresponding to the four different temperatures tested. 

What are the explanations? In summary, either polybutadiene does not 
fall into the category of materials satisfying the reduced variables law 
or some other parameter besides frequency and temperature is affect- 
ing the damping and loss modulus data. In the light of the preceding 
discussion of high dissipation results the possibility exists that 
the influence of strain softening is reflected even in the low 
(reference) dissipation results for the shear specimens. 

Since increasing frequency at constant dissipation means a reduction 
in strain, the low dissipation test data for damping is likely to 
reflect the effects of strain if these are significant. If now we 
review the effects of strain on damping for shear specimens (Fig. Ill), 
it is apparent that decreasing strain will also decrease damping. 

Thus, on the basis of strain effects alone, "low dissipation" damping 
would be expected to fall with increasing frequency. This anticipated 
trend is a plausible reason for the inconsistency of loss modulus data 
and the method of reduced variables. 

Stiffness is seen to increase with decreasing strain and, therefore, 
with increasing frequency at constant dissipation. This is a trend 
which is consistent with the method of reduced variables, one of whose 
basic tenets is that an increase in frequency has an effect correspond- 
ing to a decrease in temperature. 
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Fig. 116 FMS Deviation in Storage Modulus vs Characteristic Temp 
erature. Shear Specimen (Eight Elements) 
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Fig. 119 Method of Reduced Variables. Reduced Loss Modulus vs 

Reduced Frequency. Reference Temperature, 32 C. Char- 
acteristic Temperature, -50 C 



Before any conclusion can be drawn with respect to the validity of the 
reduced variables method for the particular material it is necessary 
to obtain further shear specimen values at amplitudes sufficiently low 
to exhibit strain independence of the material properties. In addition 
data at constant amplitude should be obtained to eliminate this influence 
of strain. 

With the present material property data for polybutadiene there is 
limited benefit in applying the reduced variables method as a method 
of combining the effects of frequency and temperature. Since data is 
available at a range of temperature values of interest it has been 
summarized and stored directly in that form. 

For application the curves of Figure 120 have proved most useful. In 
these curves storage and loss modulus at 100 and 1000 Hz are plotted 
against temperature. Values at these two frequencies may therefore be 
read for any temperature between 32 and 93 C. To calculate stiffness 
and damping at these two frequencies for a particular element the 
appropriate geometric relationships are applied, and plots of stiff- 
ness and damping against frequency are made which exhibit a straight 
line variation on log- log paper. This is how the predicted data for 
cartridge and cylindrical specimens presented in the preceding Section 
were generated. 

Consistency of Phase II And Phase III Results 

To provide a check on consistency, the material property data obtained 
from Phase II tests (as described earlier in this report) has been 
compared with the Phase III material property values of Table X for 
32 C, The comparison, presented in Table XIV shows that there is good 
consistency between the two test series for storage modulus, but that 
the loss modulus values obtained from the most recent series of tests 
are substantially higher than the Phase II values. 


TABLE XIV . COMPARISON OF PHASE II AND PHASE III VALUES 
FOR STORAGE AND LOSS MODULUS 


Test 

Series 

G' 

100 Hz 
N/m^ 

G' 

1000 Hz 
N/m^ 

g" 

100 Hz 
N/m^ 

g" 

1000 Hz 
N/m^ 

Phase II 

1.353x10^ 

2.067x10^ 

2.021x10^ 

1.684x10^ 

Phase III 

1.369x10^ 

2.189x10^ 

3.659x10^ 

2.728x10^ 


The reason for the discrepancy in loss modulus is the simplified, more 
accurate test method for material properties employed in the most recent 
series of tests. By eliminating the preload piston from the shear tests 
a source of error has been eliminated - an error which had been exaggerated 
(rather than reduced) by the "correct ions’* applied to account for it in 
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LOSS MODULUS 1000 Hz 




reducing the Phase II data*. Contrary to expectations based on an anal- 
ysis of the mechanical system, assuming rigid system elements (ref. 2), 
the effect of the piston’s presence is to reduce the apparent damping 
as described earlier in this report. Since the correction was based 
on the analysis of the mechanical system, the correction was applied 
with the wrong sign. The reasons for this anomalous effect of the 
piston are not clearly understood, but may well result from flexibility 
of such elements as the piston itself. The fact that the present test 
data for shear specimens is more reliable is not, however, in doubt. 

Comparison of Prediction And Test Data 

The test data for cartridge stiffness and damping and for compression 
element stiffness and damping at different temperatures and the com- 
parison with prediction methods, has provided a basis for assessing 
the worth, in particular applications, of the prediction methods dis- 
cussed earlier 

For cartridge stiffness and damping predictions, the available methods 
are based on static stiffness analysis. Two simple methods, the Beam- 
Column analysis and the method of Gobel predict stiffness and damping 
which generally bracket the test data. 

The difference in results of these two prediction methods is close to 
a factor of two. The limitation in state-of-the-art prediction for 
this geometry may, therefore, be summarized as follows: 

Using static analysis-based predictions, dynamic performances of a ring 
cartridge may be predicted with a band of uncertainty which represents 
a factor of two in stiffness and damping; if this geometry is to be 
employed in a particular application, it should be ensured that satis- 
factory performance can be obtained for values of stiffness and damping 
which cover this band of uncertainty in cartridge predictability. 

For cylindrical compression specimens there exists a static prediction 
method, employing a constant shape factor coefficient value O) of two. 
Based upon this method, dynamic stiffness and damping may be predicted 
using the frequency dependent material properties of storage and loss 
modulus. Alternatively, using previously generated dynamic test data 
for cylindrical cartridge specimens with various d/h ratios, a function 
has been developed empirically which provides frequency dependent values 
for shape factor coefficient. These frequency dependent coefficients 
may again be used, in conjunction with the frequency dependent storage 
and loss module, to predict stiffness and damping. 

For the Phase II test data, the static based method generates dis- 
crepancies of as much as 50 percent relative to the test data. 
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For the Phase III test data, at temperatures of 66 and 93 C, the two 
prediction methods produce discrepancies of similar magnitude. Figures 
63 through 66 show the discrepancies graphically. The worst discrep- 
ancies between the prediction methods and test data points for all 
test condition of preload and temperature are summarized in Table XV. 
For stiffness the worst discripancies vary from 8 to 25 percent; for 
damping the discrepancies are in the range 20 to 30 percent, with 
measurements lower than predicted* 


TABLE XV. CYLINDRICAL COMPRESSION SPECIMEN 

[Maximum Discrepancies Between Measured And 
Predicted Stiffness And Damping Value] 


Fig . 
No . 

Temp . 
(C) 

Pre- 

load 

a) 

Stiffness 

Damping 

ipsi 

Percent* 


Percent 

6=2 ^ 




63 

66 

2.5 

-8.8** 

+10.8 

-23.7 

-25.8 

64 

66 

5 

+12.1 

+13.4 

-28.6 

-30.8 

65 

93 

2.5 

+13.0 

+14.4 

-20.7 

-23.5 

66 

93 

5 

+22.9 

+25.0 

-21.8 

-25.2 


^Percentage referred to predicted value. 

**Negative sign indicates measured data lower than predicted. 
Positive sign indicates measured data higher than predicted . 


In terms of performance, there is nothing to choose between the two 
approaches for cylindrical compression elements. When either method 
is to be used in design of a system component, it is necessary to 
ensure that satisfactory system performance can be achieved allowing 
for uncertainties of the magnitude indicated in Table XV. 

Test Method 

The ability of the base excitation resonant mass method to provide 
stiffness and damping data for a variety of geometries and over a 
range of temperatures and dissipation levels has been convincingly 
demonstrated by the results presented in this report. It has been 
determined that the "cleanest" test condition occurs when the resonant 
mass is supported solely by the test specimen, without either of the 
lower or upper loading and unloading pistons. This simple test con- 
dition introduces no uncertainty with respect to tare stiffness and 
damping . 
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The simple test condition with no pistons can be achieved for the 
cartridge and shear specimens, without danger of putting the elastomer 
metal bond in tension. The compression specimen as presently designed 
requires the preloading effect of the piston to maintain the bond in 
compression. A modified design can be envisaged in which a self- 
preloading mechanism is employed. 

Future Directions 

The results presented in this report make clear the need for further 
work to improve fundamental understanding and applications capability. 

The influence of strain amplitude as a parameter independent of thermal 
dissipation level must be further investigated as a function of geometry 
for polybutadiene and for different materials. The present results in- 
dicate that, over a significant range of strains, stiffness and damping 
are governed by strain rather than frequency. This must be further 
tested. In addition, the conditions of strain for which the stiffness 
and damping are independent of strain must be investigated, and related 
to geometry. It is anticipated that tests at very low amplitudes as 
well as high amplitudes will be required. The means to present the 
resultant information in a form suitable for use in design must also 
be developed. The information will be largely empirical, until the 
strain softening effect becomes better understood, and test methods 
must be evolved which can rapidly define the effect for different 
materials and geometries. 

To improve understanding of the importance of thermal dissipation and 
its relationship to strain softening, thermal properties of conductivity 
and specific heat must be determined for the material being tested and 
included correctly in the thermoviscoelastic model. Direct measurements 
of temperature within the elastomer material, and comparison of these 
measurements with predictions will also provide improved understanding 
of thermal dissipation effects. 

When strain softening is better understood, and stress independent 
results are available, the method of reduced variables as a means of 
conveniently presenting elastomer properties must be further evaluated. 

As comparisons show, the uncertainty in prediction ability for a car- 
tridge and of cylindrical compression elements when static geometry 
effects are applied, while significant, may be acceptable for some 
applications. The degree of predictive uncertainty should be in- 
vestigated further for other materials. The need for, and potential 
benefits of, improved viscoelastic analysis methods must be identified 
and, if the need is shown to exist, appropriate methods of viscoelastic 
continuum analysis should be developed. 

Tests over similar parameter ranges for geometries similar to those 
presently reported should be performed using other materials, including 
materials selected for known differences in hardness and loss coef- 
ficient* These tests will provide the basis for generalization of^ or 
exception tc^ the results for polybutadiene. 
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In the case of a cartridge, the effect of a rotating load as opposed 
to a uniaxial load must be investigated. The anticipated application 
of the cartridge element is in rotating machinery where the load 
vector will be rotating. It is important to determine, under con- 
trolled conditions, whether there is a significant difference in the 
stiffness and damping under the action of rotating and uniaxial dynamic 
loads, and to establish the degree of uncertainty in prediction for a 
rotating load environment . 

An important geometry to be tested is the 0-ring. This is a compact, 
commercially-available, simple, elastomer element. 0-rings are already 
in use as seals and flexible mounts in squeeze-film dampers where their 
stiffness properties in particular are important. Their use as flexi- 
ble mounts independent of squeeze film dampers is anticipated, in 
which case their damping properties will also be important. 

When testing 0-rings, and other geometries, it is important to preserve 
the approach which has been inherent in the tests and analysis up to 
this point; that, for improved understanding, tests in shear, for 
general material properties, must be performed on the same material as 
is employed in the geometry of interest. In this way, meaningful pre- 
dictive methods which relate stiffness and damping for storage and loss 
moduli can be evolved . 

Analytical studies for representative machinery classes should be 
performed to establish the needs for damping elements and the ability 
and limitations of elastomers in filling these needs. 

Application tests should be performed to demonstrate the vibration 
control capabilities of elastomer elements. Following successful 
application tests, the questions of life and survivability should be 
investigated. 
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APPENDIX A 


THERMAL DISSIPATION IN AN ELASTOMER ELEMENT UNDER UNIAXIAL STRESS 


As a result of imposed conditions an elastomer element, in the form 
of a right circular cylinder, is undergoing time varying changes in 
temperature, internal stress and displacement. The cylinder is as 
shown in Figure 121. It is assumed that stresses are uniaxial. 

At any point 


a = a(z,t)e^^'' 

z V > / 

= U(z,t)e^'“^ 

- 1/2, t) 

where is the normal axial stress 

is the axial displacement; is temperature; 

and it is understood that only the real parts of expressions (72) and 
(73) apply. 

a and U are complex numbers containing both amplitude and phase informa- 
t ion . 

It is assumed that the amplitude functions a,U are undergoing much 
slower changes with time than and U^. 

Axial equilibrium requires: 


(72) 

(73) 

(74) 


^ (a) = 0 (75) 

where Z is the axial coordinate and the inertia of the elastomer 
material is neglected. 

The stress amplitude a is evaluated as: 

a = 3G*(Z,<o) ^ (76) 

where^G* G' (cd,T) + i (o),T); T is the cycle average temperature 
and G , G are the storage and loss moduli. T is a function of Z. 

The instantaneous energy equation is: 
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Fig. 121 Compression Loaded Cylinder 
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where S = Re (a ) Re 


BT d^T 

p C ~ = K — + S 
P dZ^ 

^ /dU 

St VdZ 


( 77 ) 


Now, assuming temperature to be changing slowly with time, we can take 
a cycle average as follows to get T: 


^.t+2rr /cD 


(78) 


thus the energy equation becomes: 

P c. 

Now, if 


St 


m , t+(2TT/0)) 

^ = K ^ \ S(t')dt 

dZ^ 2n 


(79) 


t+(2n /co) 

I = \ S(t')dt' 

^t 


(80) 


then 


, .2tt/o) 


^ = 4 




where the superscript ~ implies the complex conjugate. 
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Thus 
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the final form of the 


^ M 1 
SZ 3Z J 


energy equation is: 


p C ^ = K ^ + 

P 3Z^ 


and the equilibrium equation is: 



(82) 


(83) 


^[3G*(Z,o.)||] - 0 (84) 

It is straightforward to perform a similar development for elements 
undergoing shear deformation (as shown in Fig, 122). In this case the 
energy equation is: 


^ = 


P St 





(85) 


where U is the amplitude of shear deformation (normal to Z) and the 
equilibrium equation is: 


iZ [g* If] = 0 (86) 

To further streamline Equation (83) through (86) to a single pair of 
equations, we can write (reorganizing slightly) for shear or compres- 
sion elements: 

Energy equation: 


S_ 

SZ 



P C 


ST 

P St 


- yg 


SZ 


J 


2 


(87) 
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Fig. 122 shear Element 


193 




(88) 


where y ~ 1 ^ shear element, 3 for a compression element. 
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APPENDIX B 


DATA REDUCTION ANALYSIS 


For the Base Excitation Resonant Mass test method, test data consists 
of accelerometer and displacement sensor readings. The primary 
objective of data reduction is to interpret these signals in terms of 
stiffness and damping. The structure of this data and the method of 
reduction can best be defined in relation to a simple spring-mass- 
damper system subjected to base excitation, as shown in Figure 123. 

The equation of motion for the mass is: 


where Y is displacement of mass 

Y is displacement of base. 

We will assume that the motions of base and mass are harmonic; so that 


Y 

= Real 


(90) 

m 


L m J 


Y 

= Real 

{y* 

(91) 

0 


L 0 J 


and the equation of motion 

may be 

rewritten : 



2 „ io)t 
- CO M Y^e 


Dividing through by we obtain 


icot 


= 0 (92) 


2 * 
M Y 

(K^ + iK^) = 

o 


(93) 


where K 2 = o) C- 

Now if cp is the phase 


Y 

Y 


* 

m 

* 

o 


3 ngle between and can write 


1 Y* 1 cos cp * i 

1 Y* 1 s in cp 

(94) 

1 

• ml 


1 ^* 0 1 


(95) 


and 
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OUTPUT, y 



INPUT, y. 


Fig. 123 Model of Base Excitation Resonant 
Mass System for Data Rednction 
Analysis 
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I 

I 


(96) 


I 

I 


(K + iK ) = ~ cos CP ~ i g sin cp") 

1 2 O' cos cp — i a sin cp - 1 

|Y*| 

where a is the transmissibility Thus 


K, + iK_ = ^,- M (fit cos y i or sin cp )(cy cos CP - 1 + i g sin cd') . 

12 ^2 2 2 (^ 2 ) 

(o’ cos cp 1 ) + O' sin cp 

or 


and 


K. 


O) M ict^ a cos cp] 

^ TT 

a za cos cp + 1 


(98) 


^ M O' sin CP 

2 

O' “2a cos cp + 1 


(99) 


Thus, if we determine a, the transmissibility, or ratio of output 
amplitude to input amplitude |y* | , and tp the phase angle between 
input and output amplitudes, we can directly evaluate stiffness and 
damping and K^). 


Note that these expressions are based on the simplest one-degree-of- 
freedom model with no correction for tare damping or stiffness. This 
IS in contrast to the expressions used in Reference 2, which made 
^rrections for effects of damping in preload and load support pistons. 
The simplicity of the present treatment results from the simplicity and 
accuracy of the test procedure performed for the majority of the 
tests, in which no external restraint of any kind was applied to the 
supporting mass. Only for the compression specimens was the preload 
piston employed. These expressions are, therefore, fully accurate for 
shear and cartridge specimens and subject to some uncertainty only 
in the case of the compression specimens. 
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NCMENCLATURE 


A 


A^,A 


2 




stressed area, 


m 


2 


coefficients of power Law relationship K, = A,cd ; 

B 

K 2 = A^a-> , N/m 

coefficients of force in generalized force-displace- 
ment relationship 


B 


factor by which frequency is multiplied to get reduced 
frequency at temperature 6 , corresponding to 
temperature 6 

2 

bulk modulus, N/m 




width of shear element in direction perpendicular to 
direction of loading, m 

exponents in power law relationship = A-Cjd , = 

A 




coefficients of displacement and its time derivative 
in generalized force-displacement relationship 


C specific heat, J/kg C 

P 


C 


T1 


damping coefficient of damping element in triparameter 
model , Ns /m 


D 


diameter of cylindrical compression element, m 


da 


elemental area, 


m 


2 


df ,df 
c’ s 


elemental forces in Beam-Column analysis, parallel and 
perpendicular to Beam-Column axis, respectively, N 


D. ,D 
1 ’ o 


s 

dt ' 

de 

E 

E 


eff 


inner and outer diameters of cartridge element, m 
energy loss per cycle due to hysteresis, J 
elemental time in cycle average integral 
elemental angle 
Young*s modulus, N/m^ 

n 

effective Young's modulus for cartridge element, N/m 
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F 

f 

F, 


f 

r 

F* 

G 

G* 

g',g' 

^eff 

^eff 


K 


K 


Kti.K^2 

K, 


^2 

£,L 


force exerted by elastomt'rj N 
frequency, Hz 

maximum force in hysteresis loop, N 

factor used by Gobel to correct for shape effects in 
cartridge elements of finite length 

reduced frequency, Hz 

complex force amplitude, N 
2 

shear modulus, N/m 

complex modulus G ^ + i(/ N/m"^ 

2 

storage, loss moduli u\ sliCur. 

effective storage modulu- I- shear, N/nr^ 

effective loss modulus in shear, N/m^ 

9 

Static shear modulus ^ N/m 
reduced shear modulus, N/m^ 

height of compression element; thickness of shear 
element (strained dimension for compression and sliear 
elements ) 

(2tt /^) X cycle average of dissipation function, S, 
J/md/cyc le 

thermal conductivity- watts/cm C 

complex stiffness, 

radial stiffness of cartridge, N/m 

stiffness values in triparameter model, N/m 

stiffness, N/m 

damping, N/m 

length of cartridge element, m 
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M 

n,N 


r 




s 

T 

t 

t 


U 

U 


m 

m 




supported mass, kg 

number of elements in test specimen 
mean radius 

right-hand side of energy equation in numerical 
solution at time, t, J/m^S 

inner and outer radius of cartridge, m 

dissipation function, j/m^S 

shape factor - ratio of one loaded area to free area 
temperature, C 
t ime , s 

dummy time for cycle average integral, s 

ambient temperature in metal to which elastomer is 
bonded, C 

instantaneous temperature at axial coordinate Z, C 
axial displacement, m 

energy stored in element at maximum displacement, J 

axial displacement (instantaneous), m 

length of shear element in direction parallel to 
loading, m 

deflection of elastomer element, m 

complex amplitude of elastomer deflection, m 

absolute displacement of supported mass, m 

complex displacement amplitude of supported mass, m 

absolute displacement of base, m 

complex displacement amplitude of base, m 

maximum deflection of elastomer element, m 

axial coordinate, m 
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transmissibility , ratio of amplitude magnitudes: output 
acceleration to input acceleration 

shape factor coefficient 

shape factor coefficient for use in calculating 
stiffness 

shape factor coefficient for use in calculating 
damping 

factor equal to 1 for shear elements, 3 for compression 
elements 

time step in numerical solution of energy equation, s 

loss coefficient, 

angle in cartridge analysis 

absolute temperature in reduced variables analysis, K 
characteristic temperature for elastomer material, K 
transition temperature for elastomer 
density of elastomer, hg/m^ 
complex stress amplitude, N/m^ 

2 

normal axial stress (instantaneous), N/m 

phase angle between motions of base and supported mass 

frequency , rad/sec 
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